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Aplicação do Método Explícito-Implícito em Sistemas 
de Equações Parabólicas, não Lineares, em Domínios 
de Geometria Arbitrária 


RESUMO: 


Define-se o método explicito-implícito para a 
solução numérica de sistemas de equações às deri- 
vadas parciais parabólicas não lineares em domínios 
de geometria arbitrária e apresentam-se alguns re- 
sultados de aplicação. As vantagens relativas ao 
método dos elementos finitos e aos de diferenças 
finitas implícitas são referidas quanto aos aspectos 
memória — tempo de cálculo. 


1. INTRODUÇÃO 


A resolução de problemas em transferência de 
energia e matéria, quando formulados em termos 
das equações fundamentais, envolve tipicamente a 
resolução das equações às derivadas parciais, de 
tipo parabólico ou elíptico. No caso mais simples, o 
sistema reduz-se à equação do calor, ou da difusão: 


dg 


ot 


No caso de escoamento em meios porosos, de 
difusão com reacção química, etc., ter-se-á em geral 
um sistema de equações não lineares. Em termos de 
complexidade crescente, o caso geral é constituído 
pelas equações do movimento tridimensional de um 
fluido viscoso, em que existe transferência de ener- 
gia (calor) e de matéria. Se o escoamento é turbu- 
lento, os modelos de turbulência actualmente desen- 
volvidos para o cálculo das tensões de Reynolds, 
traduzem-se também por sistemas de equações do 
mesmo tipo, que representaremos simbolicamente 
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SUMMARY: 


The explicit-implicit method for the numerical 
solution of systems of non linear parabolic partial 
differential equations in arbitrary geometries is 
defined and some results of application are shown. 
The relative advantages in relation to the finite 
element methods and implicit finite difference 
methods are refered in what concerns memory- 
«computing time. 


O método de solução numérica do sistema para- 
bólico anterior, é também aplicável à resolução do 
problema elíptico que corresponde à estacionaridade. 


Sob o aspecto analítico, para a maioria dos 
problemas não é possível a obtenção de soluções. 
Em inúmeros casos não existe mesmo a demonstra- 
ção formal da sua existência e unicidade. 

Numericamente, vários métodos têm sido pro- 
postos, apresentando cada um o seu domínio próprio 
de vantagens, em face da natureza do problema e 
do rigor exigido para o resultado. 


Dum modo geral, podemos agrupar os métodos 
existentes em duas grandes categorias: os que 
possuem uma formulação matemática completa e a 
demonstração de estabilidade e convergência; os que 
possuem uma formulação essencialmente física, cuja 
validade é testada com resultados experimentais 
obtidos num modelo físico. A complexidade das 
equações torna muitas vezes impossível a construção 
de esquemas numéricos para os quais seja possível 
a demonstração matemática de estabilidade e con- 
vergência. Por outro lado, tal demonstração implica 
que a dimensão da malha tenda para zero, o que 
constitui um requisito impossível de satisfazer na 
maioria dos cálculos. Sucede ainda que a realização 
em computador de um dado esquema numérico 
introduz erros não tidos ainda em conta pela teoria 
actual do cálculo numérico, pelo que, a redução da 
malha pode conduzir à divergência numérica, contra- 
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riamente ao que a teoria (incompleta) permitiria 
esperar. Os esquemas baseados em formulações 
essencialmente físicas, não raro carecem na generali- 
dade que uma análise matemática mais aprofundada 
permitiria conferir-lhes, 


2. FILOSOFIA DO MÉTODO IMPLÍCITO-EXPLICITO 


O método de cuja aplicação se apresentam 


neste artigo alguns resultados, procura conjugar as 


duas vias anteriores e visa resolver deficiências 
conhecidas noutros métodos, os quais se reflectem 
entre outros nos seguintes aspectos: 


— incapacidade de tratar geometrias arbitrárias 
que podem variar no tempo (como é o caso, 
por exemplo, dos escoamentos com superficie 
livre) 

— excessivo tempo de cálculo e/ou excessiva 
exigência de memória no computador 

— facilidade de programação e adaptação a 
problemas específicos. 


A dificuldade de adaptação a geometrias arbi- 
trárias, e sobretudo a geometrias variáveis no 
tempo, é típica dos métodos implícitos e também 
dos do tipo elementos finitos. Tal facilidade existe 
nos métodos explícitos, mas as suas características 
quanto a estabilidade são deficientes. 

Em termos de memória e tempo de cálculo, o 
método de elementos finitos é em geral o mais 
exigente sem contrapartida apreciável na precisão 
mas facilidade na geometria. Os métodos implícitos 
menos exigentes em memória e tempo de cálculo 
tém geralmente grandes vantagens em problemas 
de escoamentos a números de Reynolds elevados. 
Tanto num como noutro caso, surgem com frequên- 
cia problemas de estabilidade numérica devido ao 
mau condicionamento das matrizes finais, o qual 
aumenta em geral com o refinamento da malha. 

O método implícito-explícito, aplicável a sistemas 
de segunda ordem, baseia-se numa separação da 
malha em duas submalhas, uma das quais assente 
nos nós cuja numeração é ímpar, outra naqueles 
cuja numeração é par. Nestas condições, um qual- 
quer ponto interior par fica rodeado de pontos cuja 
numeração é impar e vice-versa. O cálculo é efec- 
tuado a partir dos valores iniciais, calculando os 
valores no instante seguinte em nós pares (ou ímpa- 
res) usando um método explícito. Os valores corres- 
pondentes a esse instante nos nós impares (ou pares) 
é feito usando um esquema implícito. Todavia, se a 
equação é de segunda ordem e o esquema usa apenas 
diferenças centrais, verifica-se que usando os valores 
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calculados no passo explicito para as incógnitas que 
surgem no implícito, a matriz das incógnitas se 
reduz neste à forma diagonal, pelo que as mesmas 
se obtêm directamente, Obtidos os valores em todos 
os nós correspondentes ao instante a seguir ao ini- 
cial, o cálculo procede do modo habitual mas alter- 
nando os pontos em que se efectua o cálculo por um 
esquema explícito. Tal alternância torna o esquema 
equivalente ao de Dufort-Frankel para a equação 
de calor (que é intrinsecamente estável) embora não 
exija a memorização de dois níveis no tempo. 

Na sua aplicação prática o esquema adapta-se 
facilmente a geometria variável, é mais rápido que 
os métodos implícitos conhecidos e revela em muitos 
casos estabilidade e precisão superiores. 


3. EXEMPLO DE APLICAÇÃO 


A informação necessária para a resolução de 
um problema consiste nas condições iniciais e fron- 
teira, além das propriedades fisicas. O programa 
testa todos os pontos, identificando os que consti- 
tuem a fronteira e o interior em cada instante, 
De acordo com os resultados do cálculo, a fronteira 
pode mover-se. Na sua forma geral, o programa 
pode tratar qualquer número de equações em duas 
ou três dimensões do espaço. 

A titulo exemplificativo considera-se o caso dum 
escoamento bidimensional transiente. As equações de 
Navier-Stokes são resolvidas na forma perturbada 
descrita em [1] e [2]. 

O primeiro exemplo refere-se ao escoamento em 
torno de um triângulo, o segundo ao escoamento 


sobre uma cavidade. Em ambos os casos o valor ' 


de Re é igual a 6 e os vectores representados indi- 
cam o valor da velocidade. 
O programa existe disponível no C.T.A.M.F.U.L,. 
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Figura 1-B — Pormenor do escoamento, na 
zona posterior do triângulo 
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Figura 2-B — Escoamento na zona superior à cavidade 


Figura 2-4 — Escoamento na cavidade 
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Figura £-C — Pormenor do escoamento na zona à saída da 
cavidade 
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O método dos elementos finitos 


RESUMO 


Neste artigo faz.se uma introdução ao método 
dos elementos finitos, descrevendo as suas principais 
características que o tornam distinto dos outros mé- 
todos de solução de equações diferenciais. 


1. INTRODUÇÃO 


O método dos elementos finitos (designado aqui 
abreviadamente por MEF) foi usado pela primeira 
vez por Courant em 1943 [5], no estudo de proble- 
mas de torsão e por Prager e Synge em 1947 [14] 
no domínio da elasticidade. 

Contudo o método só se começou a tornar popu- 
lar a partir da sua redescoberta no âmbito dos pro- 
blemas de estruturas no final da década de 50 [9]. 
De facto, o método exige vasta manipulação aritmé- 
tica pelo que o seu sucesso dependeu e ainda depende 
do desenvolvimento dos computadores digitais de lar- 
ga capacidade. Não é assim de estranhar que os tra- 
balhos pioneiros de Courant e Prager e Synge não 
tenham tido de imediato a continuação que mereciam. 
Foi preciso aguardar a chegada dos primeiros gran- 
des computadores na década de 60 para que o MEF 
pudesse florescer em toda a sua potencialidade, 

Embora o método se tenha desenvolvido para 
resolução de problemas estruturais não tardou muito 
a reconhecer-se-lhe um âmbito mais geral. Assim hoje 
é considerado um método de solução de equações dife- 
rências, ombreando e competindo neste domínio com 
o método das diferenças finitas. Esta fase coincidiu 
com a «matematização» do MEF' (ver, por exemplo, 
as referências [1, 2, 3, 11, 12, 13, 15, 17] para uma 
apresentação rigorosa do método). 

O presente artigo tem um carácter de divulgação 
pelo que as preocupações didáticas sobrelevarão os 
aspectos mais matemáticos, não havendo pretensão 
alguma em se ser completo ou exaustivo. 
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ABSTRACT 


This paper presents the finite element method in 
an introductory way and describes its principal fea- 
tures which distinguishes it from other methods for 
solving differential equations, 


2. MOTIVAÇÃO 


Neste parágrafo apresentamos um exemplo mui- 
to simples para mostrar como aparece e funciona o 
MEF. Imaginemos que queríamos resolver o seguinte 
problema: calcular a deformação e as tensões numa 
viga simplesmente apoiada (ver figura 1) submetida 
a uma carga vertical uniforme p. 


Fig. 1 
Como sabemos a equação da flexão simples de 


vigas aplicável a este caso é: 


EI =p (2.1) 


em que os símbolos significam o seguinte: 


E — módulo de Young 

I — momento de inércia 

u — deslocamento vertical 

x — coordenada em relação ao apoio esquerdo 
p — carga sobre a viga 
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Uma vez conhecida a função u(x) podem cal- 
cular-se as extensões e a partir da lei de Hooke as 
respectivas tensões. 


No caso presente a integração da equação com 
az condições aos limites 
u(to) = 0, u(l) = 0 (2.2) 
não oferece quaisquer dificuldades, A solução que 


aparece em qualquer livro de Resistência de Mate- 
riais, é 


] 3 
u apa nm P. > 
24 ) 


(P—2be +xº) (2.3) 


Como é que resolveríamos este problema se em 
vez de termos condições tão simples como as que 
existem neste caso, tivessem, por exemplo, uma dis- 
tribuição de cargas mais complicada, propriedades 
elásticas variáveis (i.e. E ser função de x), geome- 
tria, variável (i.e. I ser também função de x), etc,, 
de modo que não soubéssemos integrar a equação 
diferencial (2.1)? 


A pista é-nos fornecida pela observação de que 
a solução de (2.1) é a função u(x) que minimiza o 
seguinte funcional 


| l 


: + 
J(v) = / EI Rad dx -| p(x) vdx (2.4) 
: 2 dx E 
(6) OQ 


sujeito às condições de fronteira (2.2). Quer isto dizer 
que se u for a solução de (2.1) com (2.2) como res- 
trições então 
J(u) < Jtiv) (2.5) 
para todo o função v tal que v(o)=v(l)=o e que 
seja admissível. Veremos mais adiante o significado 
preciso deste termo. Por enquanto basta nos obser- 
var que as funções v devem ser contínuas e possuir 


derivada contínua, pois caso contrário a viga não 
suportaria as cargas impostas. 


Notemos que o primeiro termo de (2.4) repre- 
senta a energia elástica da viga e o segundo termo 
o trabalho das forças aplicadas, pelo que em termos 
físicos o princípio enunciado e simbolizado em (2.5) 
mais não é que o princípio do trabalho virtual: em 
cquilíbrio o trabalho virtual das forças aplicadas é 
igual à variação de energia elástica. 


Agora temos à disposição duas vias para resolver 
o nosso problema: o da integração da equação dife- 
rencial e a minimização do funcional associado. Como 
admitimos, para efeitos de exposição, que a primeira 
não é praticável vejamos como podemos aproveitar 
a segunda. 

Uma das formas de o conseguir é propor uma 
função v(x) para o deslocamento em que haja parâ- 
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metros livres e minimizar relativamente a esses pã- 
râmetros. Concretizemos um pouco, 


Suponhamos que a função 


v(x)=a ta,x+a, x? (2.6) 
era proposta como representando aproximadamente 
o deslocamento da viga. (Notemos de passagem que 
um polinómio de grau inferior a dois não seria acei- 
tável para o caso vertente pois que, com as condições 
de apoio (2.2), à função v(x) reduzir-se-ia ao caso 
trivial v(x)=0.) 

Os valores de a, a, e a, serão determinados de 
forma a que v(x) dado por (2.6) minimize a funcio- 
nal J(v) sujeita aos constrangimentos implicados pe- 
las condições de apoio, Como 


| l 
J(v) = | Er Qa)tar— fria t ax-+a,x')dx (2.8) 


“ 


(O) O 


ou seja 


. 2 N) 
J(v) = 4EI a—p alta, -—— +a, — (2,9) 


Os constrangimentos são 


a =0, a+tal+a, E =0 (2.10) 

A pesquisa do mínimo de J(v) dado por (2.9) 
constrangido a (2.10) pode ser feita usando a técnica 
dos multiplicadores de Lagrange. Esta consiste na 
formação do funcional modificado J* (v, À) dada por 


J*(vA)=J(v)—A, (a,) — a,(a, ta, l+a, E) 2.11) 
que deve ser minimizado sem constrangimento. Os 


N's são conhecidos, como se sabe, pela designação de 
multiplicadores de Lagrange. Então devemos ter que 


“q (2.12) 


no ponto mínimo. Estas expressões fornecem o se- 
guinte sistema de equações lineares 


3 
seta, — rr = PÉ (2134) 
3 
M+A = pl (2.13b) 
| 
n= 5 (2130) 
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(2.13d) 


a+la,+Pa, — (2.13€) 

A natureza dos multiplicadores de Lagrange tor- 
na-se clara a partir da análise das equações (2.13b, c). 
A soma de À +A, é igual à carga total sobre a viga 
(pl) e lk, é igual ao momento flector relativo ao 
apoio esquerdo ou direito. (Esta ambiguidade é re- 
sultante de o problema que estamos considerando ser 
totalmente simétrico), Não é difícil pois verificar 
que os A's não são mais que as reacções de apoio. 

É vantajoso escrever as equações (2.13) sob a 
forma matricial. Teremos assim 


100 0:00 a |j=| 0 [talã) 
DR a, j 
O O SELO (-E a, pl3/3 
0 06 114 É pl 
10 00 0) 1 A pl/2 


mn E 4 = E O no JE 


Estamos na presença de um sistema de 5 equa- 
ções a 5 incógnitas: 3 «deslocamentos» a, a, e a, 
e 2 «reacções» à e A,. Para evidenciar a natureza 
distinta destas variáveis particionámos a matriz aci- 
ma em conformidade. Não é difícil ver que o sistema 
(2.14) tem solução única como aliás a situação fisica 
exige, sendo possível demonstrar que tal deve acon- 
tecer. Após a determinação dos a's temos à nossa 
disposição uma função que descreve aproximadamen- 
te (assim o esperamos) a deformação da viga. Com 
este dado podemos calcular as extensões e tensões 
que lhe estão associadas. O problema estã em prin- 
cípio resolvido. Há no entanto vários aperfeiçoamen- 
tos desejáveis. 

Em primeiro lugar não é claro qual o significado 
físico a atribuir aos coeficientes a, a, e a... Em se- 
gundo lugar as incógnitas a's e A's aparecem mistu- 
radas no sistema (2.14), Ambas estas desvantagens 
podem ser resolvidas usando outras funções em vez 
do polinómio (2.6). Consideremos que 


v=4, (xa, +y, (x)a, +y, (x)q, (2.15) 
com 
à (== E (É + (2.164) 
l l 
Ya 06) = e ( — — 1) dabiniid, 


A função v continua a ser um polinómio do 2.º 
grau com três parâmetros q,, q, e q, livres. A dife- 
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rença reside nas funções y (x) que os afectam. Con- 
centrando-nos sobre a função v verificamos que ela 


l 
se anula para x= —— ex=letomao valor um 
p 


para x = 0. Podemos dizer que ela representa um 
«modo» de deformação unitário da viga, em que os 
pontos médio e o apoio direito permanecem como 
nodos, i.e. imóveis, Similarmente para y, e y, que 
são deformações unitárias dos modo central e do 
apoio direito, respectivamente. A sobreposição destes 
três modos de deformação de acordo com (2.15) dá 
a defrmação da viga. Podemos repetir agora os cál- 
culos que fizemos atrás. Aproveitaremos contudo 
toda a simplicidade da notação e álgebra matricial 
para abreviar as expressões. Assim se representar- 


mos por (.) um vector-coluna e (. ye o seu trans- 
posto (vector-linha) e [.] uma matriz, teremos 
(2.14) escrita desta forma 


e por sua vez 


u= x) la) (2.17) 
du dy T T 
a q) = (y | (2.182) 
elo pos 
dºu d% |r r 
= |—— + (9) =(Wy') (q) (2.18b) 
dxº | dx? | 
obtendo-se facilmente por substituição, 
] 
[E = EE à 
10= [5 (Fte Va laz — 
o 
| 
— | pl Vala 2.19 
o 
Se pusermos 
l 
[K] = [EL ty”) (4) dx (2.208) 
o 
l 
(£) = | p (v)dx (2.20b) 
ba 
a expressão (2.19) escreve-se simplesmente 
l À T 
NA = a IK] ta) — tal tt) (2.21) 


A matriz [K] chamaremos matriz de rigidez e 
ao vector-coluna [f) vector das forças aplicadas equi- 
valentes. Adiante justificaremos estas designações. 
Notemos apenas que [K] é simétrica como resultado 
de (2.208). 
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As condições nos apoios escrevem-se agora 


Luc (a) =6 (ema =o (223) 


Agrupando estas duas equações podemos por que 
[B] (q) — O (2.23) 


em que [B] é no presente caso uma matriz de 2x3 


cuja primeira linha é o vector ( uloy e cuja segunda 


linha é o vector (y! $o 
O funcional modificado pela introdução dos mul- 
tiplicadores de Lagrange toma agora a forma 


J* (va) =3 (v) — (NT [B] fg) (2.24) 


em que (A) é o vector-coluna dos multiplicadores de 
Lagrange. A minimização não constrangida de 
J* (v,A) segue os mesmos trâmites que atrás, obten- 
do-se no final as seguintes equações 


[KI] (aj+[BI* [a) = [f) (2.25) 


[B] (A) —0 (2.25b) 


Estes dois sistemas podem escrever-se como 
um só, 


- —— me — —— ———— 


— ——— — — — — 


Observemos que a matriz do primeiro membro 
continua a ser uma matriz simétrica, o que não acon- 
tecia em (2,14). Isto é uma vantagem que se apre- 
ciaria devidamente quando, em vez de 5 incógnitas 
que o problema actual tem, estivéssemos a resolver 
um em que aquele número ascendia às centenas ou 
milhares. A economia de memória e de tempo de 
cálculo pode ser neste caso substancial. Após a so- 
lução de (2.26) obteríamos os deslocamentos nodais 
(q) e as reacções (A). 


Em relação à segunda desvantagem apontada 
atrás é fácil de ver que é possível rearranjar o sis- 
tema (2.26) por troca a ordem das equações de for- 
ma a poder calcular o vector (q) primeiro e segui- 
damente o vector (A). Isto significa resolver um sis- 
tema de 3 equações primeiro e depois um de 2 equa- 
ções o que é mais fácil do que resolver um sistema 
de 5 equações de uma só vez. (Outra vantagem será 
referida mais adiante). 


Fazendo o ponto da situação dizemos que esta- 
mos na posse de um método aproximado de resolver 
a equação das vigas elásticas. A maior ou menor 
aproximação dependerá do grau do polinómio utili- 
zado. Embora não tivessemos feito a demonstração 
formal deste facto, dado que qualquer função con- 
tinua pode ser aproximada com um erro arbitraria- 
mente pequeno por um polinómio de acordo com o 
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teorema de Weirstrass, aquela conclusão é plausível. 
Remetemos o leitor interessado para as referências 
citadas. 

Contudo, aqui aparece uma limitação grave ao 
método. A medida que o grau do polinómio aproxi- 
mador cresce, a matriz [K] do sistema começa a tor- 
nar-se quase singular e muito difícil de inverter, ori- 
ginando erros enormes nos valores das incógnitas (q) 
e (h), ao ponto de tornar os seus valores sem signi- 
ficado. E para isto acontecer não é preciso ir muito 
longe. Polinómios de grau 10 são suficientes já para 
denunciar este comportamento, Sendo assim o mé- 
todo de minimização da funcional usando polinómios 
como funções aproximadoras que parecia tão geral 
acaba por ser aplicável apenas a problemas de pe- 
quena dimensão. É aqui que entra em cena o MEF. 


3. O MEF UNIDIMENSIONAL 


Uma maneira de obviar o inconveniente apontado 
atrás é o de considerar a viga como sendo constituída 
por elementos cada um dos quais é dotado com uma 
função aproximadora do tipo polinomial usado em 2. 
Estas funções definidas por troços deverão possuir 
a necessária continuidade. Nascem assim os elemen- 
tos finitos. Para exemplificação suponhamos a viga 
dividida em dois elementos cada um dos quais defi- 
nido por dois nodos (ver figura 2). 


1 elem.] 2 1 elem.2 2 


] 2 3 
mm 
Fig. 2 


No presente caso há que escolher funções que nos 
nodos comuns a dois elementos sejam contínuas e 
possuam primeira derivada contínua nesses nodos. 
Os polinómios de Hermite satisfazem estes requisi- 
tos. Assim para o elemento 1 usaremos as funções 


Va (É) = 1 — 342 +283 com £=2x/1 (3.1a) 
E 
Va (É) = (é — 282 + €3) =. (3.1b) 
Ya (É) = 3 é? — 2 Fi (3.1c) 
18 
Va (é) e (ES rem 6) (3.1d) 


Teremos portanto que para este elemento 


V, = (v), ta), (3.24) 
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Notemos que o 1.º parâmetro do elemento 1 (g,), 
está associado a um deslocamento unitário do nodo 1 
e o 2º parâmetro do mesmo elemento (q,), a uma 
rotação unitária do mesmo nodo. O 3.º e 4.º parâ- 
metros estão associados a um deslocamento e rotação 
unitários, respectivamente, do nodo 2. Este caso é um 
exemplo em que há 2 parâmetros por nodo, sendo 
costume falar-se em 2 graus de liberdade por nodo. 

Para o elemento 2 teríamos uma situação similar, 
i.e., 


Vo = (v), (dq), (3.2b) 
A continuidade da função v no nodo 2, 
VV, (3.34) 
(no nodo 2) 
dv, dv, 
é (3.b) 
implicam, como é fácil ver, que 
(q), = (q), (q), = (q), (3.4) 


Estas igualdades são verificadas automaticamen- 
te se identificarmos o modo comum como um único 
nodo, Na prática isto consegue-se usando, conjunta- 
mente com a numeração dos nodos a nível de ele- 
mento, uma numeração global (ver figura 2). 

Obtivemos deste modo uma função v definida por 
elementos, ou seja, 


v=v, para u<x< 1/2 (3.5a) 


v=v, para 1/2 < x< 1 (3.5b) 

Podemos introduzi-la no funcional J(v) e repetir 
os cálculos da maneira praticada no parágrafo ante- 
rior. Surge aqui uma simplificação notável que é a 
dos integrais em (2.4) poderem ser calculados ele- 
mento a elemento e somadas as contribuições de 
cada elemento no fim, Quer dizer que obteríamos 
para os elementos 1 e 2 as matrizes de rigidez [KI], 
e [K], e os vectores forças (f), e (f),. A matriz de 
rigidez global [K]será obtida a partir das matrizes 
elementares da forma seguinte. Para os nodos 1 e 3 
apenas contribuem os elementos 1 e 2, respectivamen- 
te, enquanto que para o nodo 2 ambos os elementos 
contribuem, Esquematicamente, podemos representar 
o processo de formação da matriz global ou assem- 
blagem na forma que a figura 3 descreve. 

Obtém-se assim um sistema de 6 equações (3 
nodos X2 graus de liberdade por nodo) ao qual se 
deve juntar as 2 equações que traduzem as condições 
de apoio. É possível rearranjar as equações de forma 
a que as reacções nos apoios possam ser determinadas 
posteriormente ao cálculo dos deslocamentos (q). Isto 
permite utilizar o facto de que a matriz de rigidez 
tem uma estrutura em banda simétrica para econo- 
mizar em memória e tempo de cálculo, 
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O processo de melhoria de precisão é conseguido 
agora pelo aumento do número de elementos, em vez 
d: aumentar o grau dos polinómios usados. 


[k] (a) = [t) 


(a) (t] 


nó 2 


no3 


Fig. 3 


4. UM CASO BIDIMENSIONAL 


Neste parágrafo vamos mostrar como é que se 
pode usar o MEF quando o problema é pluridimen- 
sional. Utilizaremos para ilustração o problema da 
condução de calor em sólidos isótropos, 

Seja 92 um domínio (limitado) no plano euclidea- 
no e à a sua fronteira, que suporemos constituída 
por duas partes à, e 92, dijuntas (ver figura 4). 
A equação que rege a condução de valor é, como se 
sabe, a equação de Poisson. 


4 4 
ÇÃO seg mó (4.19) 
6x? oy* 
u=a em 39, (4.1b) 
du 
— + gu=y em à9, (4.1c) 
on 


em que os símbolos significam o seguinte: 


u — temperatura 
x,y — coordenadas cartesianas 
g — potência das fontes de calor 
a-— temperatura imposta na fronteira 029, 
n — normal exterior à fronteira 292, 
aB e y são funções de x,y mas não de u. 


As equações (4.1a) e (4.1b) traduzem as condi- 
ções de fronteira. A primeira é uma condição de 
temperatura imposta e a segunda de fluxo de calor 
imposto. 

É possível provar que a minimização do funcional 
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E ' (4.2) 
J (v) - || e + (SO) Jaa 2 f ivan 
0X , 0y ú 
o o 


- 


| Bv? ds — 2 | yYv ds 
AM, o, 


conduz à equação (4.la), desde que se tome como 
classe de funções v admissíveis a classe das funções 
contínuas que verificam a condição de fronteira 
(4.1b). Pelo facto de esta condição ter de ser imposta 
às funções v, condições de fronteira deste tipo dizem- 
-se essenciais. As condições de fronteira do tipo 
(4.1c), pelo contrário, não necessitam de ser impos- 
tas às funções v, sendo automaticamente verificadas 
pela função u que minimiza J(v). Dizem-se por isso 
naturais. 


Notemos que agora as funções aproximadoras 
só precisam de ser contínuas, podendo as suas deri- 
vadas de 1.º ordem serem descontínuas. Este é uma 
consequência do facto de que a equação de Poisson 
(4.1a) é de 2.º ordem e o funcional (4.2) envolve ape- 
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nas as 1.º derivadas. Recordemos que ne exemplo 
da viga a equação diferencial era de 4.º ordem e o 
funcional envolvia a 2.º derivada, 

A utilização de polinómios em x e y sofre das 
mesmas desvantagens que evidenciamos atrás, acres- 
cidas ainda do facto que agora a verificação de con- 
dições de fronteira do tipo (4.1b) não ser, em geral, 
nada fácil. 

Vamos ver como é que se podem construir ele- 
mentos finitos bidimensionais. O elemento bidimen- 
sional mais simples é o triângulo (ver figura 4), 
sendo o domínio £ facilmente dividido em triângulos. 
É claro que se a fronteira for curva haverá aí uma 
aproximação tanto melhor quanto mais fina for a 
malha de elementos finitos usada. No caso de uma 
fronteira poligonal a cobertura do domínio for ele- 
mentos finitos triangulares será perfeita. 

Vejamos agora como se constrói a função v(x,y) 
em cada elemento. A hipótese que se sugere como 
mais natural é o polinómio de 1.º grau 


v(x,y) =a, +ax +a,y = (p(x,y))-(a) (4.3) 


Os parâmetros (aj poderão ser calculados em 
função de valores nodais. Tomando para nodos os 
três vértices do triângulo e estabelecendo para cada 
um deles uma equação do tipo (4.3) obtém-se o sis- 
tema 


q, ni 1 X 4, 
q, 1 X, Yo 2, | (4.4) 
ds 1 x, 


X, Y, a 
em que q, representa a temperatura no nodo i. 
Em linguagem matricial 
(q) = [AJ (a) (4.5) 
Se o triângulo não for degenerado, i.e. os vér- 


tices não forem colineares, [A] é invertível e podem 
escrever 


(a) = [AJ (q) (4.6) 
Então 
v=(p(xy)) [AJA (q) (4.7) 
Pondo 
(aGey)) [AJA = (yly)) (4.8) 


obtemos uma expressão já nossa conhecida 


v(sy) = (y04y))” (a) (4.9) 
Cada um dos polinómios do 1.º grau y, (x,y) 


representa um deslocamento unitário do nodo i com 
deslocamentos nulos dos outros dois nodos. 
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É agora fácfl obter 


E E 4.10 
6X E EX (1,7) | (q) ( ) 
E EE | (x,y) V ta) (4.10b) 
0y 0v 
E conveniente por 
(94) 7 Es (x,y) — (x,y) r (4.11) 
0 
de forma que 
av ty 1? a 
Y ) +( - ) = (a) toy) toy) ta) (4.12) 


0y 


Antes de prosseguir temos que verificar se as 
funções assim construídas são contínuas quando se 
passa de um elemento a outro. Como v(x,y) é linear 
em todo o elemento, é-o também em qualquer das 
arestas. Fica aí portanto determinada pelos seus va- 
lores nos nodos (vértices) que definem essa aresta. 
Identificando os valores nodais de nós comuns a ele- 
mentos adjacentes a continuidade pretendida é asse- 
gurada. Concluimos pois estar na presença de funções 
admissíveis no funcional J(v). 

A substituição na expressão (4.2) de (4.9) e 
(4.12) permite escrever, à semelhança de (2.19), 


J(v)=(a) [om (0y)" dA (q) — 2(q) | f(y) dA 


0 
(4.13) 


* 


ta” | B(v) (WJ'ds— 2 (q)” | y (v) ds 


dn, A 


Pondo 


[K] = [com (ay) dA + faty) (y)'ds (4.148) 


Q on, 
(= [eq ass | vv as 


n in, 

obtemos para J(v) uma expressão idêntica à (2.21). 
O prosseguimento do cálculo, i.e. assemblagem, impo- 
sição das condições de fronteira essenciais, etc., se- 
gue o mesmo caminho que no parágrafo 2, pelo que 
nos dispensaremos de o repetir aqui. No final teremos 
um sistema de equações lineares cuja solução nos dá 
o valor das temperaturas nos nodos. 

Convém notar que em muitos casos os integrais 


em (4.14) não são calculados analiticamente, mas 
sim, numericamente. 


“ 


(4.14b) 
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5. PROBLEMAS MATEMÁTICOS DO MEF 


Até aqui temos vindo a expor o MEF sem preo- 
cupação de rigor matemático. Convém no entanto 
alertar o leitor para os seguintes pontos que é preciso 
ter em atenção na justificação do método: 


a) — em primeiro lugar a escolha das funções 
base (v), a sua continuidade e admissibilidade relati- 
vamente ao funcional J(v) 


b) — em segundo lugar a demonstração da exis- 
tência e unicidade da solução 


c) — em terceiro lugar a demonstração da comn- 
vergência e estimativa de limites superiores do erro 


d) — influência do uso de integração numérica, 
se for caso disso, nos resultados 


e) — influência do uso de aritmética de precisão 
finita nos resultados, 

Para a resolução destes problemas e outros de- 
vem consultar-se, por exemplo, as referências (1, 3, 
15, 17). 


6. IMPLEMENTAÇÃO EM COMPUTADOR 


É óbvio do que ficou dito nos parágrafos ante- 
riores que o MEF só é praticável dispondo de um 
programa de faça os cálculos que o método requer. 
Descreveremos em seguida a organização básica de 
um programa de elementos finitos. Assim um pro- 
grama deve permitir a: 

a) Preparação de dados. O MEF requer mesmo 
para problemas de reduzida dimensão um volume de 
dados apreciável (malha e descrição de elementos, 
coordenadas nodais, propriedades materiais, condições 
de fronteira, etc.). E conveniente dispor de rotinas 
que façam a maior parte possível deste trabalho. 

b) Verificação de dados. E indispensável uma 
verificação tão exaustiva quanto possível dos dados. 
Isto evitará o prosseguimento do programa com da- 
dos errados com os inconvenientes que facilmente se 
advinham (desperdício de tempo de computador, pos- 
sibilidade de o erro não ser detectado e os resultados 
tomados como verdadeiros, etc.). 

c) Formação da matriz e vector de forças para 
cada elemento, Esta é a primeira fase de cálculo pro- 
priamente dito. Utilizando as expressões pertinentes 
calculam-se as matrizes de rigidez e os vectores das 
forças nodais elemento a elemento. Para alguns ele- 
mentos as integrações têm de ser feitas numerica- 
mente, i.e. são substituídas por somatórios pondera- 
dos apropriados. As informações necessárias nesta 
fase e relativas às características dos elementos a 
usar devem encontrar-se incorporadas naquilo a que 
se costuma designar por biblioteca de elementos. 


d) Assemblagem. Esta é a rotina que forma o 
sistema de equações final por assemblagem das ma- 
trizes de rigidez e vectores força elemento a elemen- 
to. Duas vias podem ser seguidas. Se o problema 
cabe todo em memória a assemblagem pode ser feita 
à medida que as matrizes de rigidez e vectores de 
força forem sendo calculados. Caso a dimensão do 
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problema não permita este procedimento devem estes 
armazenar-se em bandas ou discos magnéticos para 
que a rotina de assemblagem os possa ir lá ler. 

e) Solução do sistema de equações. Um progra- 
ma de elementos finitos depende em elevado grau da 
qualidade das suas subrotinas para a solução de sis- 
tema de equações lineares. Estas devem permitir a 
solução em memória ou por blocos e devem tirar par- 
tido da estrutura (banda, simetria, etc.) da matriz 
de rigidez global. A rapidez desta fase à muito im- 
portante e a precisão dos resultados é um ponto crí- 
tico. 

f) Saídas. Após a obtenção dos valores nodais 
podem calcular-se quantidades de interesse (tensões, 
fluxos, etc.). Esse é o objectivo das rotinas de saída. 

g) 'Pós-processadores. O volume de resultados é 
geralmente muito grande e nem sempre facilmente 
absorvível na forma numérica. É assim muito útil 
dispor de saídas gráficas, em «plotters» ou «displays» 
ópticos. 

Com todos estes requisitos não admira que os 
programas de elementos finitos «a sério» atinjam as 
dezenas de milhares de instruções em FORTRAN. 

As referências [4, 6, 7, 8, 9, 10, 18] contem 
descrição detalhada de forma de aplicar o MEF a 
problemas de estruturas, condução de calor, mecáã- 
nica de fluidos, etc. e da sua implementação 
prática. 
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O método das singularidades no estudo tri-dimensional 
do escoamento em turbomáquinas axiais!!) 


SUMÁRIO 


Apresentam-se os fundamentos do método das 
singularidades para representação do escoamento em 
turbomáquinas axiais. Em particular, descreve-se a 
aplicação dos métodos da linha e da superfície susten- 
tadoras, para estudo do escoamento em cascatas 
rectilineas de pás com geometria tri-dimensional. 


1 — INTRODUÇÃO 


O estudo teórico do escoamento através das co- 
roas de pás fixas e rotativas das turbomáquinas axiais 
constitui um capítulo importante da aerodinâmica, 
devido ao vasto campo de aplicação deste tipo de 
máquinas (ventiladores, compressores, turbinas hi- 
dráulicas, de vapor e de gás). Acrescente-se que, 
por razões práticas, é em geral difícil obter resulta- 
dos experimentais suficientemente pormenorizados que 
forneçam ao investigador e ao projectista uma ima- 
gem adequada do escoamento. 

Sendo o perfil alar adoptado no fabrico da maio- 
ria das pás, não surpreende que o estudo teórico do 
escoamento neste tipo de máquinas tenha notáveis 
semelhanças com a teoria aerodinâmica da asa de 
avião, nem que muitos dos métodos teóricos desenvol- 
vidos para este caso tenham sido utilizados posterior- 
mente, com as necessárias modificações, no estudo 
daquele. 

A par de tais semelhanças, existem consideráveis 
diferenças entre os dois casos, que tornam especial- 
mente difícil o estudo do escoamento nas turbomá- 
quinas. Em particular, mencionamos a interferência 


A. F. DE O. FALCÃO () 


SUMMARY 


The use of singularity methods for the analysis 
of the flow im axial turbomachines is outlined. In 
particular, lifting-line and lifting-surface theories are 
employed to study the flow about a three-dimensional 
straight cascade of twisted blades. 


entre pás, a presença de paredes de forma mais ou 
menos complexa, e o carácter não permanente do 
escoamento devido à presença de pás fixas e móveis. 

Quando as paredes de revolução são de forma 
aproximadamente cilíndrica, um dos métodos mais 
usuais de simplificação consiste em substituir, em 
cada ponto, o escoamento real tri-dimensional por um 
escoamento bi-dimensional plano em torno dum sis- 
tema infinito de pás dispostas paralelamente umas 
em relação às outras, com perfil, passo e posição 
relativa obtidos por planificação da superfície cilín- 
drica de revolução de raio r que passa pelo ponto 
considerado (Fig. 1). Este tipo de escoamento (cas- 
cata bi-dimensional rectilínea) tem sido objecto de 
numerosos trabalhos, e existem hoje métodos teóricos 
relativamente satisfatórios para a sua determinação 
(pelo menos nos casos de escoamento permanente 
incompressível ou subsónico), além de abundante 
documentação experimental (vejam-se, por exemplo, 
[1] e [2]). 

Este método de redução do modelo tri-dimensio- 
nal a outro bi-dimensional (blade element theory) 
não tem em conta a interferência entre as diferentes 


(1) Comunicação apresentada ao Primeiro Congresso Nacional de Mecânica Teórica e Aplicada, Lisboa, Dezembro de 1974. 
(7) Professor Extraordinário do Instituto Superior Técnico, Lisboa. 
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ct, 


+ —s 


E 


Fig. 1 


superfícies de revolução (diversos valores de r), ou 
seja, parte da hipótese de que as linhas de corrente 
formam superfícies de corrente de revolução. Esta 
hipótese é em geral aceitável quando a relação r,/T. 
entre os diâmetros interior e exterior é próxima da 
unidade. O mesmo não sucede quando aquela relação 
é consideravelmente inferior à unidade, especialmente 
se a circulação do escoamento em torno de cada pá 
varia apreciavelmente ao longo desta. Neste último 
caso, o primeiro teorema de Helmholz mostra que 
existe um rasto de vórtices (trailing vortices), seme- 
lhantes aos que ocorrem no escoamento em torno 
duma asa de envergadura finita. Tais vórtices indu- 
zem um campo de velocidades caracteristicamente 
tri-dimensional, cujo estudo teórico constitui o prin- 
cipal objectivo deste trabalho. 


2 — METODO DAS SINGULARIDADES 


Consideremos o escoamento bi-dimensional pla- 
no em torno duma asa de envergadura infinita (ou 
limitada por duas paredes planas perpendiculares 
ao seu eixo longitudinal). Se a asa for substituida 
por um vórtice rectilíneo com igual circulação, a 
força de sustentação por unidade de envergadura 
não é alterada, e, a uma distância suficientemente 
grande comparada com as dimensões do perfil, os 
dois escoamentos são aproximadamente idênticos. Se 
pretendermos porém uma representação mais adequa- 
da na vizinhança da asa, teremos que recorrer a sis- 
temas mais complexos de singularidades. Para escoa- 
mento não viscoso, um método exacto consiste em 
substituir a asa por uma distribuição de vórtices per- 
pendiculares ao plano do escoamento, tendo como di- 
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rectriz o perfil da asa. Um método mais simples, que 
dá uma aproximação relativamente satisfatória para 
perfis delgados com curvatura e ângulo de ataque 
pequenos, consiste em tomar a corda do perfil como 
directriz das linhas de singularidades (vórtices, fon- 
tes e poços). Estes métodos foram posteriormente 
aplicados ao escoamento bi-dimensional em torno 
duma cascata rectilínea de pás. Em particular, citam- 
-se os trabalhos fundamentais de Schlichting [3], 
utilizando o segundo dos métodos referidos, e de Mar- 
tensen [4], baseado no primeiro, que permitem, para 
uma dada geometria, determinar o campo de veloci- 
dades. 


3 — TEORIA DA LINHA SUSTENTADORA 


Para o caso tri-dimensional, uma das simplifi- 
cações clássicas, primeiramente utilizada por Prandtl 
para o estudo da asa de envergadura finita, consiste 
em supor que os vórtices livres, que constituem o 
rasto de cada pá, têm a sua origem num único vór- 
tice ligado à pá — linha sustentadora — cuja inten- 
sidade é igual à circulação total da pá na secção 
recta correspondente. Uma simplificação adicional, 
aceitável se é pequena a deflexão do escoamento pelas 
pás (o que em geral sucede nos ventiladores e com- 
pressores), consiste em considerar os vórtices livres 
do rasto como sendo convectados pelo escoamento 
não perturbado (escoamento na ausência de pás). 
Este método de linearização do problema significa, 
para escoamento permanente (pás fixas), que os vór- 
tices livres coincidem com linhas de corrente do es- 
coamento não perturbado. No caso de pás de um 
rotor, se o escoamento relativo é permanente, os vór- 
tices livres coincidem com linhas de corrente rela- 
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tivas do escoamento não perturbado e têm formas 
aproximadamente helicoidais (Fig. 2). 


L — linha sustentadora 
VL — vórtices livres 


que, no caso geral, o perfil, e portanto a circulação, 
não são constantes ao longo das pás. As linhas sus- 


Fig. £ 


Neste trabalho, vamos considerar um modelo fi- 
sico mais simples e supor que as pás estão dispostas 
segundo uma cascata rectilínea (Fig. 3), admitindo 


Fig. 3 
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tentadoras que representam estas são vórtices recti- 
líneos perpendiculares às duas paredes planas que 
limitam o escoamento, A montante, a uma distância 
suficiente das pás, o escoamento não perturbado é 
uniforme, com velocidade U, e faz um ângulo ,/2-—-e 
com o plano que contém as linhas sustentadoras. 
Deste modo, introduzindo a hipótese de pequena per- 
turbação, os vórtices livres são semi-rectas e formam, 
para cada pá, uma folha de vórtice plana. Utilizamos 
um sistema de coordenadas cartesianas x,y,Zz (com 
vectores unitários 1,j,k) e tomamos para unidade de 
comprimento a distância entre as paredes. 

Vamos designar por U* a velocidade do escoa- 
mento e escrevemos 


U*” =U+V, (1) 


em que U é a velocidade do escoamento não pertur- 
bado e V a perturbação devida à presença das pás. 
Admitindo escoamento sem viscosidade, o termo de 
perturbação pode ser escrito com a forma [5,6] 


oo 
V =grado—i-jtan )H(y)T (z) s s(u-nt). (2) 
n=-0 


Nesta equação, 4» é uma função potencial, H(y) a 
função de Heaviside, 5(y) a função de Dirac, t o 
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passo da cascata e T(z) a intensidade da linha sus- 
tentadora (circulação em torno da pá) que se supõe 
em geral variável ao longo da pá. 

No segundo membro de (2), o termo de sumató- 
rio representa claramente a presença das singulari- 
dades. Para tornar mais óbvio o seu significado, apli- 
camos o operador rotacional. Tendo em conta que 
rot grad 4=0, resulta 


So 
rotV = ks (y)T(z) x 5 tu-nt) — 


n=- 


o 
—l(itane+j) H(y)T' (z) 3 6 (u-nt). 
n=-% 


Desta expressão facilmente se conclui que, no segun- 
do membro, o primeiro termo (proporcional a T(z) ) 
representa os vórtices ligados, e o segundo termo 
(proporcional a T'(z) ) representa os vórtices livres. 
Excepto nestas singularidades, o escoamento é irro- 
tacional. 

Para escoamento incompreensível, da aplicação 
da equação da continuidade div V = O resulta a se- 
guinte equação de Poisson para g 


oo 
lap g=secieH (y)T(z) s 8 (u-nt) — 
n=-% 


o 
—tanes (y)T(z) s 8 (u-nt). (3) 
n=- 


Entrando com as condições de fronteira d45/02=0 
para z=0, z=1, e V=grandg=0 para y-—, à 
integração analítica de (3) conduz a uma solução 
que se pode escrever com a seguinte forma (vejam-se 
em [5,6] a derivação e expressões completas da 
solução) 


Vixyo=P(z) Vxy) + V(x, y,2). (4) 


Nesta equação, V2d(x,y) é o campo de velocidades 
(apenas com componentes segundo as direcções Ox 
e Oy) do escoamento bi-dimensional que se obtém 
tomando simplesmente I[=constante=1. Conclui-se 
imediatamente que T(z) V*' (x,y) representa a solu- 
ção simplificada do problema (blade element theory) 
obtida analisando o escoamento separadamente em 
cada plano z=constante e ignorando os efeitos tri- 
“dimensionais (não nulos quando T-+constante) devi- 
dos às interferências entre os diferentes planos e à 
presença de vórtices livres. Estes efeitos tri-dimen- 
sionais estão representados pelo termo adicional 
V“ (x,y,2Z), cujas três componentes são em geral 
não nulas. 

O método clássico de utilização destes resulta- 
dos, de modo a relacionar a geometria da cascata 
com o escoamento, consiste em substituir, para cada 
plano z=constante, o campo vectorial [(z)V* (x,y) 
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— que descreve o escoamento bi-dimensional em torno 
de pás de perfil pontual e circulação T — pelo campo 
de velocidades (que designamos por Wº* (x,y; Zz) ) do 
escoamento bi-dimensional em torno duma cascata 
de pás cujo perfil e posição relativa são definidos pela 
intersecção da cascata real com o plano z=cons- 
tante. 

Ao substituir o escoamento real pela sobreposi- 
ção dos três campos de velocidades U* (x,y,2Z) = 
=U+W“(x,y;2z) + V“ (x,y,2z), há que ter em con- 
ta, quando se procura estabelecer o campo de velo- 
cidades Wº (x,y; Z) (por via analítica, numérica ou 
experimental), que este se sobrepõe a um escoamento 
não uniforme U+V*“ (x,y,Zz), e não a um escoamen- 
to uniforme como é habitual no estudo das cascatas 
bi-dimensionais. Para contornar esta dificuldade, re- 
corremos ao método utilizado por Prandtl na teoria 
da asa de envergadura finita, que consiste em substi- 
tuir, para determinação de W* (x,y;z) em cada plano 
z=constante, o campo vectorial U + V'! (x, y,z), na 
vizinhança de cada pá, por um escoamento uniforme 
fictício definido por U + vsd(z), em que vid(z) = 
= V:d(0.0,2z). Ainda na hipótese de pequena pertur- 
bação, é lícito substituir vs! (z) apenas pela sua 
componente v,º%º (z) perpendicular a U e à linha 
sustentadora (Fig. 4). 


Fig. 4 


Note-se que uma das principais dificuldades do 
problema reside no facto de que os campos vecto- 
riais Vid(x, y, Z) (ou v/Sd(z) e Wad(x, y, z) são 
interdependentes, De facto, enquanto V“ (x, y,2z) de- 
pende de T (z), o campo de velocidades W2' (x, y; Z) 
ve portanto T(z) ) depende da direcção do vector 
U+v," (z). Por outras palavras, os termos bi-di- 
mensional e tri-dimensional não podem ser calculados 
separadamente. O modo de resolução analítica e nu- 
mérica deste problema está descrito em [5,6]. No 
capítulo 5, apresentamos resultados numéricos, sob a 
forma de curvas, para cascatas com diversas configu- 
rações, juntamente com resultados obtidos por ou- 
tros métodos, para comparação. 


4 — TEORIA DA SUPERFÍCIE SUSTENTADORA 


Existem fundamentalmente dois factores limita- 
tivos da validade da teoria da linha sustentadora, tal 
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como foi apresentada acima. O primeiro refere-se à 
simplificação feita ao substituir o campo de veloci- 
dades U + Visa (x, y, z) por U + v,3! (z). Recorde-se 
que esta substituição tem em vista a utilização dos 
valores de Wº (x,y; z) obtidos pelos métodos clássi- 
cos das cascatas bi-dimensionais. O segundo factor 
está intrinsecamente ligado ao modelo físico adopta- 
do, que supõe, para cada pá, os vórtices livres do 
rasto como tendo a sua origem num único vórtice 
ligado (linha sustentadora), e não em todos os pon- 
tos da superfície da pá ou em todos os vórtices liga- 
dos que, no modelo menos simplificado da superfície 
sustentadora, representam a pá. 


Embora o método da linha sustentadora seja 
frequentemente utilizado no estudo tri-dimensional do 
escoamento em turbomáquinas axiais, a ordem de 
grandeza dos erros acima apontados, inerentes ao 
método, não é bem conhecida. Para esclarecer tal 
problema, iremos recorrer a um modelo físico mais 
perfeito e complexo, que consiste em representar cada 
pá por uma superfície de vórtice a ela ligada (super- 
fície sustentadora) e não por único vórtice. 

Neste trabalho, iremos supor que se trata de pás 
suficientemente delgadas para que se possa conside- 
rar a sua espessura como aproximadamente nula. 
Aliás, admitindo pequena perturbação, é aceitável 
analisar separadamente e sobrepor os efeitos devidos 
à espessura das pás e os que estão relacionados com 
a força de sustentação (ângulo de ataque, curvatu- 
ra). Este procedimento é clássico na teoria dos perfis 
aerodinâmicos. 


O método da superfície sustentadora consiste em 
substituir a pá por uma superficie — superfície de 
esqueleto — em geral curva, definida pelo lugar geo- 
métrico das linhas de esqueleto dos perfis das suas 
secções rectas, e na qual estão situadas as linhas 
de vórtice representativas da pá. Continuando a admi- 
tir pequena perturbação, substituimos aquela superf- 
cie por um plano — plano de base — tanto quanto 
possível próximo dela, e no qual existem os vórtices 
ligados à pá e também os vórtices livres do rasto 
(Fig. 5). A relação entre o campo de velocidades, 
induzido pelos vórtices, e a geometria da cascata tra- 
duz-se pela condição de que a componente da veloci- 
dade normal à superficie de esqueleto é nula. Iremos 
admitir que a superfície de esqueleto e o plano de 
base são suficientemente próximos para que se pos- 
sam substituir os ângulos do escoamento calculados 
num ponto daquela superfície pelos correspondentes 
valores calculados na projecção do ponto sobre o pla- 
no de base, 

Na análise do problema, consideramos apenas 
pás de corda constante e usamos um sistema de coor- 
denadas cartesianas x,y,Z, com o eixo Oz coinci- 
dente com o bordo de ataque duma pá (ou com 
a sua projecção sobre o respectivo plano de base). 
Os planos de base das pás ficam assim definidos 
por x—-y tan e=nt (n=0+1,+2....). Para es- 
tudar o escoamento em torno da pá de referência 
n = 0, é conveniente introduzir duas outras coordena- 
das i;=xsinetycoseen=-—-xcose+gyeine, 


TECNICA 437 


Fig. 5 


que, juntamente com z, formam outro sistema de 
coordenadas cartesianas. Para a pá n = 0, os bor- 
dos de ataque e de fuga (ou as suas projecções sobre 
o plano de base) ficam definidos por £= 0 e E =c, 
n = 0, sendo c a corda da pá, 

Seja então q (£,z; E) a velocidade, induzida num 
ponto P (£,7=0,z) do plano de base n=0, por uma 
cascata de linhas sustentadoras de intensidade 
P(z, :') cuja localização é definida pelas coordenadas 
t=t' e m=nt cose (n=0,+1,+2,...). Utilizando 
resultados referidos no capítulo 3, teremos 


qts 2, E)= VI(£-£)cose, (E-£) sin É z], 


em que V(x,y,z) é dado por (4). 

Passando agora ao modelo físico que nos inte- 
ressa, consideremos uma cascata de superfícies sus- 
tentadoras, ou seja, uma distribuição continua de li- 
nhas sustentadoras em cada plano de base, com in- 
tcnsidade T(z,t) por unidade de comprimento ao 
longo da corda. A velocidade total Q(£, z) induzida 
por todos os vórtices (ligados e livres) no ponto 
P(£ n=0,2) é então obtida por integração 


Qu =fa(gzg) ar. (5) 
0 


Seja n=f(£,z) a equação da superfície de esque- 
leto da pá de referência n=0. Admitindo pequena 
perturbação (Q<U), a condição de tangência do es- 
coamento à superfície de esqueleto pode escrever-se 
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(6) 


em que Q, é a componente de Q na direcção perpen- 


dicular ao plano de base, e q é o ângulo (que se su- 
põe pequeno) entre a direcção do escoamento não 
perturbado (vector U) e o plano de base. A resolução 
analítica e numérica deste problema foi tratada em 
[5, 6). 


5 — RESULTADOS NUMÉRICOS 


Nas Figs. 6, 7 e 8 apresentam-se resultados nu- 
méricos obtidos com base no modelo da linha susten- 
tadora (abreviadamente designado por LS) e no mo- 
delo da superfície sustentadora (SS). Para completar 
a comparação, as Figs. 6 e 7 mostram também re- 
sultados da teoria bi-dimensional (2D) (blade element 
theory) e de mais dois modelos que utilizam simpli- 
ficações frequentemente adoptadas no estudo de tur- 
bomáquinas. 


pra e SS SE 7 
—— superfície sustentadora 
1.8- ——-- linha sustentadora e. 
«- placa deflectora AV Sid 
——-— plano deflector / Pá 


— = — bi- dimensional 
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O primeiro destes dois modelos, que designamos 
por plano deflector (PLN), resulta duma simplifica- 
ção do modelo da linha sustentadora em que se admi- 
t2 que o passo t é suficientemente pequeno (em ri- 
gor infinitesimal) para que o termo tri-dimensional 
de velocidade (designado no capítulo 3 por V34) se 
possa considerar como invariável com a coordena- 
da x. Esta hipótese simplificativa corresponde, na 
geometria cilíndrica, ao conhecido modelo axi-simé- 
trico do disco deflector (actuator disc) (veja-se por 
exemplo [7]). 

O segundo modelo, que denominamos placa de- 
flectora (PLC), pode-se considerar como resultando 
do modelo da superfície sustentadora por uma sim- 
plificação semelhante à que acima referimos para 
PLN. Ao contrário de PLN, PLC tem em conta que 
a deflexão se processa ao longo duma distância fini- 
ta, entre o plano que contém os bordos de ataque e 
o plano que contém os bordos de fuga das pás (ou 
seja entre 0O<y<t cos g). Este modelo, quando trans- 
posto para geometria cilíndrica, é análogo ao que <e 
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designa na literatura de língua inglesa por «through 
flow analysis» (veja-se por exemplo [7]). 

Os resultados numéricos apresentados foram cal- 
culados para pás de compressor, de corda constante 
e perfil parabólico (flecha relativa 0,05). A fim de 
evidenciar os efeitos tri-dimensionais, foram esco- 
lhidas pás torsas, com ângulo de torsão variando 
linearmente ao longo do comprimento da pá e valor 
total (de extremidade a extremidade) igual a 10º. 
Para direcção do escoamento não perturbado foi 
tomada a tangente à linha de esqueleto da pá junto 
ao bordo de ataque da secção z=1/2 (equidistante 
das extremidades). As Figs. 6 e 7 mostram resulta- 
dos obtidos com diversos valores de t (passo), c 
(corda) e r (ângulo de disposição das pás). Recor- 
de-se que o comprimento das pás é unitário. 

A Fig. 6 representa a variação do coeficiente de 
sustentação C, (z) = 2T, t2/Uc (T,., — circulação 
total em torno da pá) e a Fig. 7 mostra a varia- 
ção do coeficiente adimensional o. = 20 E (6 z)/U 


(Q (£, z) =perturbação de velocidade à superfície da 
pá). Note-se que C, é de facto um coeficiente de 
pressão linearizado. Finalmente, na Fig. 8, apresen- 
ta-se a distribuição da componente Q, (£,z) na su- 
perfície da pá, dada por SS. Na Fig. 6, as escalas 
foram escolhidas de modo a que as rectas represen- 
tativas dos resultados dados pela teoria bi-dimensio- 
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nal tenham inclinações iguais, Deste modo torna-se 
mais evidente a influência da geometria da cascata 
sobre o comportamento tri-dimensional do escoa- 
mento. 

A Fig. 6 mostra que, com c e e constantes, os 
efeitos tri-dimensionais (em valor relativo), dados 
por SS, não são apreciavelmente afectados pela varia- 
ção de t. Tornam-se porém mais importantes com o 
aumento de c, para e e t/c constantes (Fig. 6c), e 
com o aumento de e para t e c constantes. Em 
todos os casos, verifica-se que os efeitos tri-dimen- 
sionais (ou seja o desvio em relação aos resultados 
dados por 2D) previstos por SS (o mais realista dos 
modelos considerados) são maiores do que os previs- 
tos por LS, PLC e PLN. Destes três modelos, PLN 
é o que dá pior aproximação, como seria de esperar 
(veja-se também Fig. 7). Mantendo t/c e e constan- 
tes, a Fig. 6c mostra que a variação de c não afec- 
ta os resultados de 2D e PLN, e ainda que LS se 
aproxima de SS à medida que c diminui (como su- 
cede na asa de avião. 

Tem interesse analisar os resultados de LS jun- 
tamente com os de PLC, a fim de comparar os erros 
introduzidos na avaliação dos efeitos tri-dimensionais 
quando a simplificação (a partir de SS) consiste em 
tomar a corda como nula (caso de LS) ou o escoa- 
mento como invariável com a coordenada x (caso 
de PLC). A Fig. 6 mostra que LS melhora, em rela- 


e i 
fa N | 
/ face convexa 


7 ESTO Ae SS Dr, 
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ção a PLC, à medida que t aumenta e c e e dimi- 
nuem. Verifica-se também que, na vizinhança das 
paredes z=0 e z=1, as curvas dadas por LS têm 
uma inclinação mais correcta (tomando SS como 
padrão) do que as obtidas com PLC, Por exemplo, 
no caso em que t=c=0,5, e=45º (Fig. 6c), LS dá 
melhor aproximação do que PLC perto das paredes, 
sucedendo o inverso na zona central da pá. 


6 — CONCLUSÕES 


Verificou-se que a teoria bi-dimensional (blade 
element theory) é inadequada para descrever escoa- 
mentos em que a geometria das pás é claramente 
tri-dimensional. Para estudar os efeitos tri-dimensio- 
nais, foram utilizados quatro modelos, sendo o da 
superfície sustentadora o mais perfeito. A precisão 
do: resultados dos outros três modelos varia conside- 
ravelmente com a geometria do sistema, sendo o pla- 
no deflector, de todos, o que dá resultados menos 
satisfatórios. 

Note-se que, embora os resultados apresentados 
se refiram a escoamento incompressível, é possível 
estendê-los a escoamento subsónico utilizando a co- 
nhecida transformação de Góthert [8]. 


Este trabalho foi realizado no âmbito do Projecto 


de Investigação TLE-3 do Núcleo de Estudos de En- 
genharia Mecânica do Instituto de Alta Cultura, 
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Escoamento potencial bidimensional numa cascata de perfis 
rectos — método da transformação conforme 


RESUMO 


Apresenta-se a análise do escoamento potencial 
bidimensional numa cascata inclinada de perfis rectos 
a partir da transformação conforme do escoamento 
em torno do círculo. Obtêm-se as expressões da velo- 
cidade, ângulo de ataque e forças que actuam sobre 
03 perfis. Apresentam-se e discutem-se os resultados 
numéricos obtidos para a distribuição de pressão na 
região da cascata, ângulo de deflexão, coeficiente de 
sustentação e coeficiente de interferência. 


1 — INTRODUÇÃO 


O escoamento potencial bidimensional em casca- 
tas de pás assume particular interesse, numa intro- 
dução ao vasto campo dos escoamentos em turbomá- 
quinas axiais, pois é susceptível de poder ser estudado 
teoricamente através da utilização de diversos méto- 
dos analíticos. 


A situação do problema bidimensional no caso 
geral do estudo tridimensional do escoamento em 
turbomáquinas axiais e as hipóteses simplificativas 
que a ele conduzem, sendo devidamente analisados 
em [1], transcendem no entanto o âmbito deste tra- 
balho. 


Os métodos analíticos que se utilizam no estudo 
bidimensional do escoamento em cascatas de pás po- 
dem ser divididos em dois grupos principais: 


a) Método das singularidades, que consiste em 
substituir cada pá por uma distribuição de fontes, 
poços ou vórtices, o que conduz, em geral, a uma 
equação integral cuja resolução implica a utilização 
de métodos numéricos. 


b) Métodos que envolvem tranformações confor- 
mes e que permitem passar do escoamento situado 
num plano ((£mn) para o escoamento numa cascata 
de perfis situada num outro plano z(x,y), utilizando 
ou não planos intermédios. De notar que este tipo de 
métodos, porque resolve analiticamente de uma for- 
ma completa o problema, fornece a possibilidade de 
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A. 1. DE ALMEIDA SARMENTO 
J. A. FALCÃO DE CAMPOS 


SUMMARY 


The theoretical analysis of the two-dimensional 
potential flow about a staggered cascade of straight 
blades is presented based on the conformal transfor- 
mation of the flow about a circle. Expressions are 
derived for the velocity, angle of attack and forces 
acting on the aerofoils. Numerical results for the 
pressure distribution about the cascade, deflection 
angle, lift coefficient and interference coefficient are 
also presented and discussed. 


testar a precisão dos métodos numéricos empregues 
na resolução pelo método das singularidades. 

Um dos métodos compreendidos neste segundo 
grupo consiste em transformar o escoamento em tor- 
no do círculo no escoamento numa cascata inclinada 
de perfis rectos sem espessura e foi pela primeira 
formulado por Weinig [2], encontrando-se tratado 
com relativo detalhe em [1] a [5]. 

Da formulação que a seguir se apresenta é possí- 
vel extrair alguns resultados numéricos com interesse, 
em particular a distribuição de pressão na região da 
cascata e sobre cada pá, para vários valores da razão 
passo-corda, inclinação da cascata e ângulo de ata- 
que. De salientar ainda que este método constitui a 
base para a determinação de coeficientes de interfe- 
rência para perfis em cascata, que encontram larga 
utilização no projecto de ventiladores e compressores 
axiais (veja-se por exemplo [6]). 


2 — FORMULAÇÃO ANALÍTICA 


Consideremos a transformação de Joukowski: 


(1) 


e ti)=t+ 


2 
id 
” 


O potencial complexo z representa o escoamento sem 
circulação em torno do círculo de raio a, e resulta 
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da sobreposição de um escoamento uniforme z = « 
com um escoamento produzido por um dipolo situado 
no centro do círculo. 


Se considerarmos que o escoamento uniforme é 
equivalente ao escoamento produzido por um dipolo 
situado no infinito podemos, por separação e localiza- 
ção adequada da fonte e do poço que o constituem, 
determinar um escoamento em torno do circulo de 
raio a, susceptível de se transformar no escoamento 
numa cascata de perfis rectos sem espessura. 


Desloquemos assim, a fonte e o poço de que se 
compõe o dipolo no infinito para os pontos |; = — a/k 
e (= a/kK respectivamente, sendo k um número real 
positivo inferior a 1, A fim de manter o círculo como 
linha de corrente do escoamento torna-se necessário 


É 
x" -— ' 
e vorLIce 


vortice 
E 


Para obter o caso geral de uma cascata inclinada 
é necessário introduzir no plano do círculo dois vór- 
tices do mesmo sentido nos pontos — a/kK e ak e dois 
vórtices de sentido contrário nos pontos — ak e a/k. 
Virá para o potencial complexo: 


I | 7 Ç + ak , É je ak 
Za (e ih) In talk Dad +(g + ih) In Ea HE | 


— ak 
(3) 


sendo h a intensidade de cada vórtice, Verifica-se 
que o transformado do circulo é um conjunto infinito 
d: segmentos de recta (perfis) desfazados de g se- 
gundo o eixo imaginário e de h segundo o eixo real 
(fig. 1). 


vortice 


Fig. 1 


separar o dipolo situado no centro do circulo deslo- 
cando a fonte para o ponto $;=— ak e o poço para 
o ponto ( — ak, Obtém-se assim para o potencial 


t—+ajk 


- age ak 


E +ak - 
a | (2) 


L— ak 


sendo g o caudal associado a cada uma das singula- 
ridades. 


Interpretando z como a função que define uma 
transformação conforme do plano ; no plano z, é 
evidente que o circulo se transforma numa série infi- 
nita de segmentos de recta paralelos ao eixo real 
e distanciados de g correspondendo aos infinitos valo- 


res da função logaritmica afastados entre si de 2qi. 
Sendo E um número real é evidente que a cascata 
assim obtida não possui inclinação, 
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Introduzindo a distância d (passo) entre bordos 
d: ataque e o ângulo de inclinação 5, temos h = 
— d.siná e g = d.cosg donde, desprezando a constante 


i(m/2—8) 


aditiva d/2.e + podemos escrever 


8 | ido a+ kr e £ + ak 
E E CR SS RES 1 
n + e 1% ak | (1) 


2 a — ki 

O parâmetro k estã relacionado com a razão 
passo-corda d/c. Com efeito, a localização dos bordos 
de ataque e de fuga de um perfil obtém-se da condi- 
ção dz/dj = 0, ou seja 


-1B[ k | “| | 
e ç2 pin | - 
a—ko ak 


= 


fazendo (=a 6 obtém-se 


a-+ k 
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E 
| + k? 


tgo = tgb . (5) 


Designando por O, e 6, as raízes desta equação 
tem-se O, =0,+r e fazendo + —a e 


io 
=ae ” podemos obter a distância c =x, —x, en- 
tre os bordos de ataque e de fuga (pontos correspon- 
dentes a $, e &,), ou seja a corda. Virá, depois de 


um cálculo relativamente extenso 


Vkt + 2k?cos2j + 1 + 2kcosB ê 


e - | cost In 
DR o DEP A “ 
d 


T |] — k? 
4 9ksins 
pap |. (9 
Vki + 2kcos28 + 1. 


O escoamento no plano z (plano da cascata), 
para que origine forças sobre os perfis, deverá ter 
ângulo de ataque não nulo e circulação diferente de 
zero ao longo de qualquer contorno que englobe no 
seu interior um perfil. 

O problema consiste, pois, em determinar no 
plano ; um escoamento com circulação em torno do 
círculo que se transforme por (4) no escoamento 
atrás referido. 

É óbvio que o escoamento em questão deverá ser 
um escoamento potencial conservando o círculo como 
linha de corrente e possuindo singularidades apenas 
nos pontos (=+a/k e (=+ ak, O escoamento defi- 
nido por 


É = Vem e n 


o 


de” cms 8 ak” i3 |] 
+ n 
a— k? t— ak 


(7) 


satisfaz estas últimas condições mas dá origem no 
plano z a um escoamento uniforme de velocidade V 
com ângulo de ataque nulo. 

Consideremos então o escoamento definido por 


vd[ —i(B+ta), a+k 
F = a e a” 
qe | (8) 
|— ak 


No plano z somos conduzidos a um escoamento sem 
circulação com ângulo de ataque q. Note-se que este 
escoamento não é deflectido pela cascata. Com efeito 
a velocidade no infinito a montante é 


-— dF 


e, atendendo a que o ponto correspondente a — o no 
plano | é (=—-a/K temos 
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Analogamente obtém-se para a velocidade no infi- 
nito a juzante 


Witdroo) dF º di ia 
CO | == E = o arauto 
dd lj- ak: dz t=ap=Ve 


ou seja, V(-00)=V( +00). 


Importa pois sobrepor ao escoamento (8) um es- 
coamento com circulação, Logicamente este pode ser 
obtido colocando dois vórtices do mesmo sentido em 
t=ta/k e os correspondentes vórtices de sentido 
contrário em (4 = ak. Virá para o potencial com- 
plexo deste escoamento 


ikK 12 — a2k? 

F=— Inh———— 9 
sendo K a intensidade de cada vórtice. Evidentemente 
a circulação em torno do círculo será "T=2K. Sobre- 
pondo os dois escoamentos o potencial será F + F, 


e a circulação determina-se fazendo intervir a condi- 
ção de Joukowski nos bordos de fuga dos perfis: 


d(F + F,) 
CE liap (10) 


” ” 


o que conduz a 


9ksina 


norma e o (11) 
Vk4 - :k2cos95 1 
Se determinarmos para o escoamento final as 
velocidades no infinito antes e depois da cascata 
achamos 


| iz SksinaV i'r/2— B) 
Wi ÇÃO | ——esaemeees É , 
Vk4s + 2k?coszp + 1 
(12) 
i 9ksinaV — i(x/2 o 
Víto)=Ve' + — + id 
Vks + vk?cos2p+1 


Decompondo segundo as direcções paralela e normal 
à frente da cascata virá 
QksinaV 


V (c0)=Vsinte+BD)+-——]J][][[ - 
t yV ki + 2k?coszp+1 


ec 2ksinaV 
V. (pão) = Vinte E |) =——————e—essmetri , 
t Vk4 + zk?cos2p + 1 


V (o)=V (,00)=Veos(2+8) (13) 
n 


n 
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O mesmo resultado poder-se-á obter aplicando a equa- 
ção da continuidade a ambos os lados da cascata e 
calculando a circulação da velocidade para um con- 
torno constituído por duas rectas paralelas aos perfis 
distanciadas de g e fechado por dois segmentos de 
recta paralelos à frente da cascata situados no infi- 
nito a montante e a juzante. 
O ângulo de ataque a' (ângulo da velocidade em 
= — 9), tendo sido alterado, passa a assumir o valor 


2k 


cd sina +- cosfsina 


cosa + R sinfsina 


e o ângulo do escoamento deflectido tem o valor: 


R ek ii 
— | Sina — cosísina | 
a sin 09 
cosz — —— sinfsina 
R B 


com R = Vka + Zk?cos2p + 1 


Finalmente a aplicação do teorema da quanti- 
dade de movimento conduz às componentes X e Y 
da força exercida sobre um perfil: 


k 
X = — 4;V2d —— sin? | 
AE cais 
(16) 
k 
Y == 45V2d —— sinacosa 
R 
Apresentam-se a seguir alguns casos particulares 
com interesse: 1.º A = 0 (cascata sem inclinação). 


Neste caso a equação (5) pode ser resolvida em 
ordem a k obtendo-se 


Cc 9k Cc <= 
k-==tanh | — À — = — 
an a á ) Como R tanh ( aa ) 


as componentes da sustentação tomam os valores 


c ” 
NX = — 2Vid tanh |— — 2 
an e 7) sima, 


+) 15 
Y =2.Vid tanh |— — | gi 
an ( ia ) sinicosz 


Se c/d é suficientemente pequeno será tanh 7 Ea 
2 
Cc 
- TUE e obtêm-se as conhecidas expressões das 
componentes da sustentação para um perfil isolado: 


oV2 
A — Ir — E sinêz 
2 


(18) 


V2 
? pV : 
Y “28 "a € sinzcosz . 
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oo TR ONES (14) 


Se, pelo contrário, c/d »1 então tanh Es e) =1 
e teremos: 

XxX oV2 

— 6 É “— “fiiiita , 

d po) 

(19) 

É pit, 

— 4 sinzcosz ” 

d 2 


Ainda no caso 8 = 0, interessa considerar o que se 
passa com a deflexão introduzida pela cascata no 
escoamento. Os ângulos a' e q” são dados por 


A Cc = 
tg =| 1+tanh Fr =) |rez , 


tga” t= tah (5 5) 
ga” = — tan a a )|'8 » 


o que acarreta para a deflexão 


tg. (20) 


—. 


Se c/d > 1 então tga” = 0, isto é, a veloci- 
dade no infinito a juzante da cascata é aproximada- 
mente paralela aos perfis, sendo o ângulo de deflexão 
igual ao ângulo de ataque. Podemos considerar, de 
acordo com (20), que, no que respeita a deflexão 
produzida no escoamento pela cascata, o factor 


(1 —e Cm/ q constitui o que se pode chamar a efi- 


ciência da cascata. 

2.º B = 90º (perfis alinhados ao longo do eixo 
real). 

De (5) obtém-se 


a E o (21) 


3.0 c/d-—> O (perfil isolado). 


k +. 
Neste caso tem-se — -» k e k — Ed  obten- 
R d 4 


do-se para a sustentação 


sV2 
X=— 9 “—c sin? : 
2 
(22) 
EV , 
Y=2z a c sinacosa . 
No que respeita aos ângulos q' e q” verifica-se 


que «,a” -—> a O que corresponde ao caso de deflexão 
nula, como seria de esperar. 
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3 — RESULTADOS NUMÉRICOS 
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Nas fig. 2, 3, 4 apresentam-se, para várias geome- 
trias de cascatas e para diversos ângulos de ataque q”, 
as linhas de corrente do escoamento e a distribuição 

P— P.oo 
po 1/22 V 2.00 
na região da cascata de perfis. De notar que o escoa- 
mento quase não é perturbado a uma distância da 


de pressão adimensionalizada E 


ordem de grandeza da corda, quer a montante quer 
a juzante da cascata, e que as linhas de corrente são 
já praticamente paralelas na zona do bordo de fuga. 
Como seria de esperar todas as linhas de pressão 
constante encontram os perfis entre o bordo de ata- 
que e o ponto de estagnação. 


Ú 


No que se diz respeito à localização do ponto 
de estagnação as fig. 5, 6, 7 mostram que: 


— Não é grandemente afectada pelo ângulo de 
inclinação à da cascata; 

—'Tem mais influência o valor do módulo do 
ângulo de ataque a' que o facto de ser positivo ou 
negativo; 

— É fortemente afectado pelo valor da relação 
passo-corda d/c. 

É também de referir que o bordo de ataque é um 
ponto de velocidade infinita, pelo que o valor da pres- 
são tenderá para — o e o bordo de fuga é um ponto 
de velocidade finita, em geral não nula. 

Nas fig. 8, 9 apresenta-se a variação do ângulo 


—e— us um o e e e 


] 
] 
] 
4 
| 
| 
| 
| 
| 


“- 
- 


A E 7 


Fig. 5— Cp sobre o perfil para pários valores de 8 a'-15º e d/c=1 


o EE oa... 


de. Gu gue 


Fig. 6— Cp sobre o perfil para vários valores de d/c w'=15º e B=45º 
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Fig. 7— Cp sobre o perfil para vários valores de a', 8=45º e d/c=1 
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Fig. 8— Am cao —a” em função de «' para vários valores de B e dic -1 


90 
RR 3 
+ a) Ng 
o « 
- És [+ 
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a 
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Fig. 9— Am =a'—-a” em função de a' para vários valores de d/c e B=40º 
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de deflexão áq=a-a” com ângulo de ataque q 
para diversas geometrias de cascata. Delas se con- 
clui que para uma dada cascata, quando o ângulo 
de ataque se aproxima do ângulo de inclinação ou do 
seu complementar, o ângulo de deflexão cai rapida- 
mente. Fora dessa zona a relação entre os dois ân- 
gulos é aproximadamente linear. 

A fig. 8 mostra que o ângulo de inclinação não 


O [e q 
“ 0 «150 - ] 
AU 2140 maço a e 


—— mo mm o .— e 


Fig, 10 — 
9. C, 


o o o 
170 - 9) - f 0) i 
qem. que Gr” ue e que que qe qe emo Ba 


Nas fig. 10 e 11 apresenta-se o coeficiente de 
sustentação C, em função do ângulo de ataque para 
valores de d/c e do ângulo de inclinação, O coeficien- 


te de sustentação é C = em que L é a 
L 


1 2;V?2 
— -— — 


força de sustentação eV " =1/2 (V Sd cê Je ave- 


em função de «', para vários valores de 58 e d/c=1 


E | dm 


Fig. 11— € r em função de a', para vários valores de d/c e 8=40º 


tem influência apreciável sobre a inclinação das cur- 
vas na zona aproximadamente linear, afectando po- 
rém os limites dessa zona. Na fig. 9 vê-se, como é 
ds esperar, que para valores de d/c pequenos o ân- 
gulo de deflexão é aproximadamente igual ao ângulo 
de ataque. Contudo para valores de d/c grandes, é 
muito mais pequeno. 
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locidade média. A fig. 10 mostra que, especialmente 
para valores elevados do ângulo de inclinação A, 
o coeficiente de sustentação C, cai mais rapidamen- 
to com o ângulo de ataque para valores positivos 


deste do que para valores negativos. 
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A fig. 12 representa a variação do coeficiente 


Cc 


asc, 
de interferência K, = sendo OC 


Lcasc. 
Lisol. 


relativo ao perfil em cascata e Chisol relativo ao 


perfil isolado, com d/c para vários valores do ângulo 
de inclinação da cascata. 

De notar que, no caso limite de d/c = 1 e 8=90º, 
a cascata transforma-se num perfil de corda infinita, 
pelo que o seu coeficiente de sustentação, e portanto 
Ko serão infinitos. O coeficiente de interferência é 
utilizado no projecto de turbimáquinas axiais para 
corrigir o coeficiente de sustentação obtido para um 
perfil isolado, sobretudo no caso de pás com pequena 
curvatura (Ventiladores e compressores axiais). 
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em função de d/c, para vários valores de 8 
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Programa para resolução duma rede de actividades 
pelo método do caminho crítico 


SUMÁRIO 


Embora existam já programas para resolução 
dum diagrama CPM-tempo por computador, a sua 
utilização além de onerosa não era didáctica. Por 
estas razões foi desenvolvido no CTAMFUL, para 
apoio da Divisão de Termodinâmica Aplicada, um 
programa para a resolução deste tipo de problemas. 


A correcta elaboração e planificação dum pro- 
jecto e sua execução, torna necessário o conhecimento 
completo da sequência e interligação entre todas as 
partes componentes. Com a evolução da técnica e 
consequente aumento da complexidade dos projectos, 
os antigos métodos de coordenação das diferentes 
actividades, tais como o gráfico de Gantt, tornaram- 
-se impotentes para realizar um perfeito controle 
do projecto nas suas diversas fases. 

Esta situação apresentou-se primeiramente nos 
Estados Unidos da América quando do projecto PO- 
LARIS, onde a marinha dos U. S. A, necessitava dum 
controle dos diferentes contratos para o programa 
de construção de mísseis e quando da construção pela 
Nemours Company dum conjunto de fábricas de pro- 
dutos químicos, No primeiro caso foi criado um grupo 
de trabalho encarregado de elaborar um sistema ca- 
paz de avaliar matematicamente as estimativas dos 
diferentes empreiteiros e determinar a data provável 
do término do contrato. Este novo sistema passou a 
ser denominado por Técnica de Avaliação e Revisão 
do Programas ou PERT (Program Evaluation and 
Review Technique). No segundo caso a Rand Corpo- 
ration lançou o método de Planeamento e Programa- 
ção de Projecto ou PPS (Project Planning and Sche- 
duling). 

Estava-se no período de 1956-58. Posteriormente 
ambos estes sistemas foram desenvolvidos, o segundo 
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SUMMARY 


None of the existing programs for solving CPM- 
-time problems is suitable for teaching purposes, to 
the author's best knowledge as they are all available 
in machine language only. Therefore, a new CPM- 
-time solving program was developed at CTAMFUL, 
with teaching purposes in mind, and has already been 
used by the Applied Thermodynamics Department 
students. 


dos quais no Método do Caminho Crítico ou CPM 
(Critical Path Method). 

No caso de um projecto, para uma análise correc- 
ta dos problemas, torna-se necessário saber o tempo 
que as diversas etapas constituintes levarão e as res- 
pectivas interdependências. Assim, se o custo não 
for um factor importante, deve-se atender à forma 
de proceder que nos leve a um tempo mínimo de con- 
clusão da obra em estudo. Deve-se pois atender espe- 
cialmente ao encadeamento de actividades que deter- 
minam o tempo máximo de conclusão da obra. Esta 
sucessão de actividades é o caminho crítico pois qual- 
quer atraso no fim de uma actividade da cadeia vai- 
-Ssº repercutir no atraso da data mais cedo do evento 
final ou seja num acréscimo de tempo da conclusão 
da obra. Desta forma de todas as actividades que 
são necessárias estarem concluídas para se poder ini- 
ciar nova ou novas actividades, deve-se atender 
essencialmente àquela que termina mais tarde pois 
é esta que irá condicionar o tempo fim das activida- 
des sucessoras e sobretudo no caso desta actividade 
fazer parte da cadeia crítica, pois ela irá influenciar 
directamente o tempo final de conclusão da obra. 
As outras actividades cujo tempo de execução pre- 
visto faz com que a sua conclusão se dê antes de sur- 
gir a possibilidade de se passar a uma fase seguinte 
ficam com uma certa margem de tempo de reserva. 
Nestes casos a actividade poder-se-á prolongar para 
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além do tempo previsto, sem ir influenciar o tempo 
total da obra, desde que a referida margem não seja 
excedida, Esta, para uma dada actividade é função 
da conclusão da antecessora e convém portanto saber 
entre que valores se situa o tempo de reserva disponi- 
vel, já que o problema a ter sempre em conta é que 
uma dada actividade não pode terminar para além 
dum tempo que provoque um aumento do tempo mais 
cedo do evento fim do sistema de actividades, Por 
outras palavras pode-se dizer que, como cada evento 
tem duas datas (data mais cedo e data mais tarde) 
e que uma actividade tem o seu início num evento 
e o seu fim noutro, temos quatro tipos de margens 
ou folgas que são a diferença entre as datas possíveis 
de início e fim, menos a duração da própria activida- 
de, (Ver a figura 1). 

Para a resolução destes problemas, foi desenvol- 
vido no CTAMFTUL um programa em FORTRAN uti- 
lizando o sistema DOS. Este programa calcula, dado 
um sistema qualquer de actividades, todas as cadeias, 
as cadeias críticas, os tempos mais cedo e mais tarde 
d: cada evento e respectiva folga, as margens total, 
livre, dependente e independente e traça a rede do 
sistema de actividades, ou apenas algumas destas 
questões, consoante o interesse do utilizador. 


Folga Civre 
Foiga folal 
Folga ince fen ctenhe 


Folga ale ben afeade 


No caso do exemplo aqui apresentado o sistema, 
de actividades era o representado na figura 2. 


Fig. & 


Foram pedidos, como se pode ver pelos resulta- 
dos obtidos e apresentados na figura 3, todas as ca- 
deias de actividades, as cadeias críticas, as margens 
e o traçado gráfico da rede do sistema de actividades. 

Neste sistema, dado o número reduzido de acti- 
vidades, foi possível na execução do traçado gráfico 
da rede, desenhar todas as ligações entre eventos. 
No caso de sistemas em que o elevado número de 
actividades torna impossível o traçado de todas elas, 
o programa cria conecções numéricas entre os even- 
tos, como se pode observar na fig. 4. 
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RESUMO 


A necessidade de representar de modo cómodo o 
enorme volume de dados gerados na solução numé- 
rica das equações às derivadas parciais que regem o 
escoamento de fluidos e o elevado custo do equipa- 
mento existente no mercado levou ao projecto e 
construção de um «display» óptico e ao desenvolvi- 
mento do correspondente «software» de operação. 
E este o trabalho que se descreve no artigo. O equi- 
pamento realizado, centrado sobre problemas de apli- 
cação, além do seu baixo custo, permitiu evidenciar 
as amplas possibilidades que tal «display» oferece 
em inúmeros outros domínios. 


|-—- MOTIVAÇÃO E HISTÓRIA DO PROJECTO 


1. Com o desenvolvimento dos métodos de cál- 
culo numérico no C.T.AJM.F.U.L. para a resolução de 
equações às derivadas parciais não lineares que re- 
gem o escoamento e a transferência de calor num 
fluido em domínios bi e tridimensionais, cedo se reve- 
laram dificuldades na interpretação do enorme vo- 
lume dos dados numéricos a que um cálculo corrente 
pode conduzir, sobretudo tratando-se de processos não 
estacionários em que interessa estudar pormenores 
que à partida se desconhece como e quando ocorrem. 
Efectivamente, mesmo no caso relativamente simples 
d> um escoamento bidimensional isotérmico descrito 
por uma malha de 20x20 pontos, para cada instante 
no tempo existem 2x400 valores para as componen- 
tes de velocidade e 400 valores para a pressão. Im- 
primir todos os valores a cada instante, além de im- 
plicar um custo proibitivo de papel, torna natural- 
mente difícil a interpretação. Por outro lado, em fase 
de estudo de novos métodos ou situações, podem ocor- 
rer anomalias ou instabilidades que se impõe detectar 
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SUMMARY 


The need for display representation processing 
numerical data generated on the numerical solution 
of partial differential equations which govern the flow 
of viscous fluids, and the very expensive equipment 
available on the market for visual display motivated 
the development of the prototype whose hardware 
and software are described. 

The work done stressed the importance of 
project development centered on applications and 
showed the cheapness and wide range of applica- 
bility of the equipment built. 


e corrigir. Gera-se assim a situação, algo paradoxal, 

2 um avanço na capacidade de solução de problemas 
complexos tornar mais difícil a sua interpretação, 
pelo volume de dados. O paradoxo é, afinal, apenas 
fruto duma inadequada representação espacio-tempo- 
ral das soluções numéricas. O método óbvio é con- 
densar tais soluções em representações gráficas tão 
próximas quanto possível, na sua visualização, do 
fenómeno físico em estudo, Num fluido, em experiên- 
cia de laboratório, é prática corrente visualizar o es- 
coamento introduzindo partículas no seu interior e 
seguir a sua trajectória. O fenómeno é depois quanti- 
ficado por um número limitado de medidas. 

A solução numérica tem, sobre uma experiên- 
cia de laboratório, em inúmeros casos, a enorme 
vantagem de ser mais rápida, permitir uma fácil 
e ampla variação de parâmetros característicos, 
de ser mais económica e, sobretudo, de ter em conta 
variáveis típicas do fenómeno real que é impossível 
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reproduzir em experiência de laboratório devido a 
factores de escala. 

Fazendo suceder no tempo tais representações 
gráficas, obter-se-á um filme (tipo desenho animado) 
do fenómeno, cuja análise se torna desse modo seme- 
lhante à de uma experiência de laboratório, ou em 
escala real, com a vantagem de ser em inúmeros 
casos a única possível. 

Como é óbvio, a representação gráfica animada 
«dilui» a quantificação fina do fenómeno em benefi- 
cio da apreciação global e da acentuação das carac- 
terísticas salientes. Todavia, uma vez obtida a carac- 
terização global, a quantificação fina obtém-se ime- 
diatamente fazendo imprimir os valores numéricos 
correspondentes às variáveis requeridas nos pontos 
e instantes desejados. 


2. O caso anterior é meramente exemplificativo. 
Em física e engenharia as possibilidades que um 
traçador de gráficos oferece são inúmeras consti- 
tuindo um instrumento altamente valioso de investi- 
gação e projecto. Um «display» óptico oferece, em 
relação ao traçador corrente de gráficos em papel, 
a vantagem da velocidade (que permite a animação 
directa) e do mais baixo custo. 

As enormes vantagens que um novo instrumento 
oferece num domínio ou em múltiplos domínios de 
trabalho não são, todavia, razões bastantes para que 
se disponha de recursos para o adquirir. Tal foi o 
caso do €C.CU.L. e C.T.AMF.U.L. Porém, a ligação 
estreita entre os dois, a enorme vantagem que a exis- 
tência de tal instrumento representa e a crença ge- 
neralizada entre nós (na altura em que o projecto 
se iniciou) de que a construção de tal equipamento 
ultrapassava largamente as possibilidades nacionais 
com os meios humanos e materiais disponíveis, cons- 
tituiu desafio bastante para quem tinha dessas capa- 
cidades uma visão menos derrotista, para atacar o 
problema. Os parâmetros em que o mesmo tinha de 
ser abordado e resolvido estavam, à partida, clara- 
mente delimitados: 


a) O sistema deveria ser compatível com o com- 
putador IBM 360-44 do CC.U.L. mas off- 
“line, devido às restrições impostas pelos con- 
tratos de manutenção que na altura não 
interessava ultrapassar. 

b) A visualização deveria ser feita num oscilos- 
cópio de que o C.T.A.M.F.U.L. dispunha. 

c) Os dados deveriam ser gravados pelo compu- 
tador em banda magnética standard. Esta 
seria posteriormente transportada para um 
desenrolador de bandas que accionaria 
«display». 

d) A traçagem de gráficos, a variação de in- 
tensidade luminosa e o accionamento de uma 
máquina de filmar deveriam ser suportados 
por um «software» de tal modo concebido que 
um utilizador corrente do CXC.U.L. o pudesse 
utilizar sem necessitar qualquer preparação 
ou conhecimento especial. 
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O ter-se previsto o accionamento duma máquina 
de filmar destinava-se apenas a futuros utilizadores 
sem acesso directo ao «display», ou para a elaboração 
de filmes didácticos. Efectivamente, dispondo do »dis- 
play» e da banda gravada, a visualização no oscilos- 
cópio é imediata e pode repetir-se quantas as vezes 
desejadas. 


A realização dos objectivos anteriores exigia um 
desenrolador de bandas... que não havia recursos para 
comprar. Foram, todavia, obtidos de sucata compo- 
nentes de unidade IBM desmanteladas. 

A tarefa inicial, de que se ocuparam Pedro G. de 
Oliveira, A, Neves de Carvalho e colaboradores do 
C.T.A.M.F.U.L. foi a decifração dos circuitos e a 
reconstrução de uma unidade de bandas. A partir 
dela, foram projectadas e realizadas as interfaces. 
A filosofia seguida no projecto é descrita em II. 

Da concepção e desenvolvimento do «software» 
necessário ocupou-se M. Quirino com a colaboração 
do programador Legateaux Martins. Este trabalho 
é descrito em III. 

Logo que o conjunto ficou operacional, o mesmo 
passou a testes de utilização. O primeiro, e mais na- 
tural, foi um exemplo concreto de Mecânica dos Flui- 
dos, o qual começara por motivar todo o trabalho 
dada a importância de que se revestia. A primeira 
aplicação esteve a cargo de J. J. Amarante dos San- 
tos e é descrita em IV. 


O projecto, na sua realização, conheceu naturais 
vicissitudes e levou alguns meses a concretizar-se. 
O custo final foi largamente inferior ao de uma uni- 
dade adquirida além do que possibilitou em aquisição 
de conhecimentos, tecnologia e... à vontade, que de 
outro modo nunca teriam sido obtidos. 


Realizado o protótipo, muitas outras possibilida- 
des foram encaradas e serão ou não desenvolvidas 
consoante os eventuais utentes. Refira-se, nomeada- 
mente, que o «display» pode ser ligado on-line ao com- 
putador IBM 360, dispensando assim a passagem pelo 
suporte intermédio das bandas magnéticas. 


Entretanto, o aparecimento no mercado de unida- 
des de bandas relativamente baratas permite encarar 
a construção e comercialização de um «display» a 
preços que se presumem altamente competitivos. 
Todavia, essa alternativa ultrapassa o âmbito do 
C.T.A.M.F.U.L. tanto mais que a equipa de «software- 
-hardware» responsável pela concretização do «dis- 
play» se encontra hoje integrada no Departamento 
de Informática da Universidade Nova de Lisboa por 
cujo arranque foi responsável o autor destas linhas. 
Finalmente, importa assinalar o objectivo de aplica- 
ção concreta que presidiu à realização a qual se re- 
flecte na concepção particular do protótipo. Tal orien- 
tação teve consequências extremamente benéficas. 

Não se reclama para o protótipo nenhuma par- 
ticular sofisticação técnica, Reclama-se, sim, a vas- 
tíssima gama de problemas que o mesmo permite 
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tratar, utilizando apenas os meios estritamente neces- 
sários em claro detrimento da exibição em que são 
férteis tantas realizações neste domínio. 


dd D. -D. 


ll — HARDWARE 


1 — Dados base do projecto 


A unidade construída é do tipo normalmente 
designada por «display X,Y gerador de vectores». 
As suas características principais são: 


— Funcionamento compatível com qualquer tubo 
de raios catódicos de osciloscópio (tamanho 
usual do ecran 8x10 cm). 

— Precisão de 256 pontos equidistantes por eixo 
(para um eixo de 10cm o valor da menor 


1 
divisão de escala é de —— mm). 
2,56 


— Geração por «hardware» do segmento de recta 
compreendido entre dois pontos consecutivos 
correspondentes ao valor de dois pares de 
coordenadas fornecidas pela unidade de ban- 
das. 

— Controle do feixe luminoso em três intensida- 
des: máxima, média e feixe apagado. 

— Sinal permitindo o controle por «software» do 
disparo de uma máquina de filmar. 

— Número máximo de pontos distintos por ima- 
gem sem causar problemas de persistência 
óptica: 1000 pontos. 


A unidade é constituída pelo desenrolador de ban- 


BLOCO 2 GAP BYTES 


DE CONTROLE 


actividade de diálogo com a equipa de «harware» fol 
altamente positiva. 


2 — Descrição do modo de gravação das bandas 
magnéticas do computador IBM 360/44 


Como o suporte de informação é a banda magné- 
tica, a filosofia da interface foi condicionada pela 
forma como a informação nela está gravada e pela 
velocidade de desenrolamento da unidade de bandas. 
Para melhor situar o problema descrevem-se sucin- 
tamente as características técnicas das unidades exis- 
tentes no CCUL nomeadamente: velocidade, organi- 
zação da informação em banda e modos de gravação 
possíveis. 

A velocidade de desenrolamento da unidade de 
bandas é de 37,5 polegadas/segundo. 

A informação é gravada em 9 pistas, havendo 
uma cabeça de leitura de nove canais que permite o 
acesso a nove bits em paralelo, sendo 8 bits de infor- 
mação e o nono de controlo de paridade. O controlo 
de paridade é ímpar (em cada conjunto de 9 bits que 
aparecem em paralelo tem de haver sempre um 
número ímpar de «1» lógicos). 

A organização «hardware» do computador é tal 
que o número mínimo de bits a que se tem acesso 
é 9, A cada conjunto de 9 bits, 8 de informação mais 
um de paridade, chama-se byte. 

Os bytes são escritos em banda em blocos de 
informação separados por espaços de fita, chamados 
GAP's, sem qualquer conteúdo informativo. O tama- 
nho dos blocos de informação é variável sendo a sua 
dimensão controlada por «software». De modo mais 
preciso, depois do último byte de cada bloco existem 
dois bytes de teste da informação do bloco. Só a se- 
guir é que aparece o «gap» propriamente dito. 

Esquematicamente a organização é 
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das, pelo osciloscópio funcionando no modo X,Y e 
pela interface. 

Paralelamente com o projecto de «hardware» de- 
correu o projecto de «software» de investigação em 
computação gráfica que realizou todo o «software» 
necessário à utilização cómoda da unidade e cuja 
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3 — Filosofia dos dois modos possíveis de gravação 
em banda 


O princípio base de registo da informação em 
bandas magnéticas é a variação de fluxo magnético. 
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As unidades IBM têm dois sistemas principais de 
registo de informações designados por phase encoder 
(codificação por fase; (FP, E)), e NRZ&I (no return 
to zero inverter). 

O primeiro destes sistemas permite uma densi- 
dade de gravação de 1600 bytes por polegada (1600 
bpi) e a forma de registo é tal que um «1» lógico 
corresponde a uma variação positiva do fluxo magné- 
tico e um «o» lógico corresponde a uma variação 
negativa. 

Quando da leitura, a tensão induzida aos termi- 
nais da cabeça de leitura tem a forma de uma sinu- 
soide modulada em frequência em que cada máximo 
de tensão corresponde a um «l+ lógico e cada minimo 
corresponde a um «o» lógico. Atendendo à densidade 
de gravação e velocidade de desenrolamento deduz-se 
que a frequência do fluxo informativo é de 60 KHZ. 

Exemplifiquemos com a seguinte sequência de 
informação 00110101001. 


d Bl | O) 


FLUXO MAGNÉTICO o 


TENSAO 


— E RS SS A O TO O TOO O OO o Oo 


O outro sistema admite uma densidade de grava- 
cão de 800 bpi sendo a frequência do fluxo informa- 
tivo de 30 KHZ. 


Neste aistema, o registo da informação lógica 
£0» não corresponde a nenhuma variação de fluxo 
magnético e o registo de «l» lógico corresponde à 
variação de fluxo positiva ou negativa. Na leitura, 
a um máximo ou um minimo de tensão, corresponde 
sempre a «l> lógico; os zeros lógicos correspondem 
a uma tensão nula, 


Com fins comparativos desenha-se os diagramas 
da escrita e leitura em NRAI do mesmo conjunto de 
informação, 
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4 — Filosofia geral do projecto 


A primeira opção c. 
modo de gravação, 

Dado que em phase encoder se torna muito difi- 
cil conseguir sincronizar o sistema devido não só à 
alta densidade de gravação como também ao facto 
do sinal gravado só possuir informação útil em ins- 
tantes determinados dependentes do sinal de sincro- 
nismo geral do computador, cedo se concluiu que este 
sistema de gravação não era susceptível de ser adop- 
tado. Pelo contrário, em NRZI sempre que há um 
máximo ou um mínimo do sinal lido estamos em pre- 
sença da codificação de um «lil» lógico. Este facto, 
conjugado com a paridade impar, permite distinguir 
os intervalos de tempo em que existe informação 
útil na banda, pois nesses instantes pelo menos uma 
das pistas terá a codificação de um «1» lógico. Quan- 
do todas as pistas contiverem zeros lógicos num dado 


ser feita foi a escolha do 


ESCRITA 


LEITURA 


instante pode-se concluir que esse não é o momento 
de leitura correspondente à informação válida, 

Depois de fixado o modo de gravação, e portanto 
a velocidade de fluxo informativo, tentou-se optimi- 
zar os valores dos seguintes parâmetros: 


— Máxima precisão do gráfico; isto é: refina- 
namento máximo do número de pontos por 
eixo, distinguíveis pela vista humana, tendo 
como dado o tamanho usual de um ecran de 
osciloscópio (8x10 em). 

— Máximo número de pontos distintos possíveis 
de representar sem que surjam problemas de 
variação de persistência, 
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FLUXO 
MAGNETICO 


TENSÃO 


A interdependência entre eles é manifesta se 
atendermos a que o aumento de sensibilidade por eixo 
implica um maior número de bits para codificar a in- 
formação e como a frequência é constante (30 009 
bytes/segundo) o número de pontos possíveis de re- 
presentar por imagem sem «flicker» diminui. 

Numa tentativa de optimizar estes parâmetros 
chegou-se à seguinte solução: 

A informação correspondente a um ponto do 
ecran é codificada em três bytes consecutivos sendo 
o 1.º (8 bits) para o eixo dos X; o 2.º (8 bits) para 
o eixo dos Y e o 3.º (8 bits) destinado a controlar 
a luminosidade do ponto e o disparo de uma máquina 
fotográfica. 

Deste modo o número máximo de pontos por eixo 
é de 2º = 256. Através do byte de controle é possível 
regular a intensidade luminosa do ponto por actuação 
no feixe de electrões do tubo podendo-se distinguir 
três intensidades: feixe apagado, intensidade média 
e» intensidade máxima. 

Por serem utilizados 3 bytes por ponto a fre- 
quência de geração de pontos distintos no ecran é de 
10 KHZ. Por «hardware» é traçado o segmento de 
recta compreendido entre dois pontos consecutivos 
(X, Yy c,) (X,, Yo C,.) sendo a intensidade luminosa 
do segmento controlada pela informação contida no 
byte de controlo c.. 


5 — Descrição do Dispositivo 
Para efeitos de sistematização da descrição o 


dispositivo pode ser reduzido ao diagrama de blocos 
a seguir apresentado. Nem todos os blocos mencio- 


ESCRITA 


— .— — .— uu o o — — 


LEITURA 


o — — "—— — "— — —. .— — — — —————— — o —— o 


nados são estruturalmente independentes mas do pon- 
to de vista lógico podem-se considerar como tal. 


Os nove sinais analógicos (correspondentes a 1 
byte de informação) presentes à saída das nove cabe- 
cas de leitura da unidade de bandas são primeiro 
tratados numa malha passiva a fim de adaptar o 
cabo de transmissão e de compensar em frequência 
as cabeças. Depois de compensados, cada um dos nove 
sinais é amplificado e depois digitalizado. Após este 
tratamento o fluxo de informação desdobra-se em 
duas vias. Por um lado é feita uma disjunção lógica 
dos nove sinais, por outro lado os oito bits de infor- 
mação são armazenados em flip flops do tipo RS. 

Esta memorização intermédia tem por fim evitar 
os inconvenientes resultantes do desalinhamento das 
cabeças de leitura. Assim que aparecer à entrada 
de um dos biestáveis um bit de informação ele é 
armazenado podendo ser depois processado sob o co- 
mando do sinal de sincronismo do sistema. Deste 
modo a informação a partir destes biestáveis estará 
sincronizada. Resta analisar como é obtido o sinal de 
sincronismo. O flanco de subida do primeiro sinal 
que aparecer à saída do circuito de disjunção faz dis- 
parar um monoestável de largura ajustável. A lar- 
gura do sinal de sincronismo assim gerado determina 
a tolerância que é admitida para as cabeças. Quanto 
maior for o tempo deste monoestável maior será a 
imunidade do sistema ao desalinhamento das cabeças. 
O limite superior deste tempo terá que ser tal que 
não possa ser lida nova informação sem o monoestá- 
vel estar de novo na sua posição estável, No flanco 
de descida do sinal assim gerado é lida a informação 
contida nos biestáveis RS, fazendo-se, em seguida 
o «reset» dos mesmos. O sinal de sincronismo actua 


TÉCNICA 437 


014095071950 


79n4131 
sv5asdv> 5r 


DV SW Ido 


sw 10 


Disidniãao DO SILAS 


mIDVEITA 2100 500 OvIvnimiTta 


ELA 


aq 


HiwnNtS 500 


OvIvsas 


W)isvaDDDs vNDDTA VO DHVvASia 


SIND IN] O YeNES 


Ag 5SOSTNdA] 


DwJvm a: 


voloorwny 
TWLtota 
OYSBIANOS 


So 


534408 


A 00 
O YYnIs avIvi 
ed é sogros VS19ONVNV v31907vny “INOSINIS 30 Sivnis 
PET YLISIa 3a ss 
Iv IIHOmI i Sivnv3 & Sivnvo 6 0 
: OVSUIANOI HI44NE SLig 500 309 


Ov)v Z17vs lota 


OVÍVZINONINTS 


1 99 
QvIvr 
“INOHINIS 
30 
H344N8 


valsoTenNy 
WuIsta 
QYSHIANOI 


v590Teny 
ovdvzlsoma 


vw Ivai sida 


277 


TECNICA 437 


ainda o sistema multiplexor que permite comandar 
nos instantes correctos as memórias tampão corres- 
pondentes ao X, Y e €C. Este sinal é também tratado 
de modo a detectar-se o início de «gap» com o fim 
de inibir os comandos das memórias tampão e repo- 
secionar o sistema de multipexagem de modo a que 
o próximo comando volte a ser para a memória tam- 
pão correspondente ao eixo dos X. 


O comando designado por «clock do c» permite 
a passagem simultânea da informação correspondente 
ao terno ordenado para a entrada dos conversores 
digitais analógicos onde a informação digital é con- 
vertida numa tensão analógica proporcional à infor- 
mação binária. Após a interpolação analógica, os três 
sinais resultantes vão finalmente atacar os amplifi- 
cadores horizontal e vertical e o circuito de controle 
de intensidade do feixe do osciloscópio. 

A descrição detalhada do dispositivo do ponto de 
vista electrónico e lógico é feita no relatório final 
que se encontra em vias de publicação, 


HW — SOFTWARE : PGBD — «PACKAGE» GRÁFICO BA- 
SICO PARA DISPLAY 


1 — Breve nota histórica 


No processo histórico de desenvolvimento de pe- 
riféricos de saída gráfica a solução dum «graphical- 
-display» comandado «off-line» por uma banda magné- 
tica representa o sucessor do «plotter» a que se acres- 
centa a capacidade de modificação da imagem e 
consequentemente de movimento da mesma. 

Enquanto que para fazer a análise dum processo 
dinâmico utilizando a saída por «plotter» se era obri- 
gado a representar a variável tempo como uma di- 
mensão espacial, o que se traduzia no desenrolar con- 
tínuo do papel para a representação do fenómeno ao 
longo dum determinado intervalo de tempo ou na 
representação apenas de instantâneos do processo em 
momentos privilegiados previstos pelo programa, 
aqui, a possibilidade de modificação da imagem ao 
longo do tempo de simulação do processo permite ao 
utilizador «assistir» ao desenrolar do fenómeno numa 
escala de tempo por ele escolhida. 


Actualmente, mesmo em centros dispondo de 
meios de computação gráfica mais sofisticados e com 
capacidade de interacção, mais de 80% da utilização 
desse equipamento é feita em saídas deste tipo, o que 
significa que as reais capacidades desses sistemas 
não estão sendo exploradas quer porque os problemas 
tratados não necessitam de facto delas e são tipica- 
mente problemas com saída gráfica apenas, quer 
porque a maior parte dos programas de aplicação 
que são utilizados provêm de adaptações de progra- 
mas inicialmente concebidos para «plotter» e utili- 
zação em «batch», 

Rosenfeld, [5], classifica a Computação Gráfica 
de acordo com o modo de representação dessa mesma 
informação à entrada e à saída do sistema. Admite 
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para tal que a informação gráfica é representada ou 
por uma imagem ou por uma descrição que tanto 
pode ser um conjunto de números (caso das coorde- 
nadas dos pontos duma curva) como uma fórmula 
matemática (algoritmo) que descreve a forma de a 
obter. 


ita Imagem Descrição 
- k 
= Tratamento Geração de 
É de imagens imagens 
= 
= Análise «Tudo o 
f de imagens resto» 
A 


Fig. 1 — Classificação de Rosenfeld 


As aplicações a que o «graphical-display off-line» 
di solução situam-se segundo esta classificação na 
categoria de problemas que usam da Computação 
Gráfica o subconjunto relativo à geração de imagens 
apenas. 

Note-se que esta classificação não é sempre uti- 
lizada e outros autores tais como Newmann e Sproull 
[6] e Prince, [7] e nós próprios em [8] apenas con- 
sideram a Computação Gráfica a partir da introdu- 
ção da interação, i.e. com diálogo Homem/Máquina, 


2 — Concepção dum Package Gráfico Básico para 
Display 


Quando se tem um programa de aplicação que 
gera imagens para serem apresentadas num ecran 
como meio de saída (on-ou off-line) existe uma etapa 
intermédia que é a da constituição da «display-file> 
a qual contém a informação gráfica a sair na forma 
d: instruções que comandam o movimento do feixe 
luminoso sobre o ecran, 


ECRAN 


PROGRAMA 
DE 
APLICACAO 


DISPLAY- 
FILE 


Para que o utilizador dum domínio de aplicação 
possa estar liberto dos pormenores técnicos associa- 
dos a essa etapa deve poder dispor duma linguagem 
de programação que contenha uma extensão gráfica, 
l.e., que permita ao utilizador o envio de imagens 
para o ecran através de instruções do mesmo tipo 
das outras saídas. 
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O package gráfico base para «display», (PGBD) 
aqui apresentado é um conjunto de subrotinas que 
permitem ao utilizador FORTRAN obter um número 
qualquer de imagens, tendo como únicas limitações, 
o número de pontos de cada imagem não exceder 1000 
e o tempo total de representação não exceder o tempo 
da leitura continua da banda magnética pela unidade 
de leitura. 

As imagens são construídas pelo utilizador na 
forma de um conjunto de pontos que descrevem um 
caminho que o ponto luminoso deverá seguir, e a in- 
tensidade luminosa que o feixe deve apresentar entre 
cada dois pontos. 

Toda a construção de ficheiros de imagens assim 
como o seu tratamento e escrita na banda magnética 
em forma compatível com o hardware» serão do 
ponto de vista do utilizador realizadas automatica- 
mente e as áreas de memória de trabalho reservadas 
sem a sua intervenção. 

O terminal gráfico a que se destina este package 
é constituído por um leitor continuo de banda magné- 
tica, uma unidade de controle e um ecran de raios 
catódicos. 


QO]— |» | —O) 


Cada instrução da «display-file» é um triplo de 
valores (3 bytes) que representam respectivamente as 
coordenadas (x, y) do ponto a que se referem e 
atributos da imagem associados a esse ponto. 

Uma. instrução faz o feixe deslocar-se da posição 
onde se encontra para uma nova posição sobre o 
ecran, podendo o trajecto ser percorrido pelo feixe 
com 3 possíveis intensidades luminosas (apagado, in- 
tensidade média e máxima). 

O terceiro byte contém ainda informação que 
permite decidir sobre o disparo ou não de uma má- 
quina fotográfica assim como outras informações de 
controlo do «hardwares». 

O PGBD assume que a origem do referencial 
estã colocada no canto inferior esquerdo do ecran e 
a unidade de controle assume que a informação cor- 
respondente às coordenadas X e Y pertence ao inter- 
valo (0,255). 

Dado que o G. D. não dispõe de dispositivo de 
refrescamento, isto é, não repete automaticamente a 
imagem representada tantas vezes quantas as neces- 
sárias para que o fósforo do ecran retenha a imagem 
por um tempo suficiente para que a vista humana a 
construa, o refrescamento é feito pelo «software» 
segundo a indicação que o utilizador dá para o tempo 
da permanência da imagem. 


3 — Descrição do PGBD 


O conjunto de subrotinas que o PGBD põe à dis- 
posição do utilizador pode classificar-se em dois tipos. 
Um que se ocupa da organização da «display-file» 
e o outro da construção das imagens. O primeiro 
é constituído por instruções do tipo OPEN e CLOSE 
cujo funcionamento pouco interessa ao utilizador ex- 
cepto em que é o veículo de valorização dos atributos 
das imagens através dos parâmetros, O segundo gru- 
po interessa mais directamente ao utilizador pois cons- 
titue a via de descrição das saídas gráficas do pro- 
blema (construção de imagens). São apenas duas as 
subrotinas deste tipo, PLOT e LINE, a primeira mo- 
vimentando o feixe luminoso apenas de um ponto 
para outro directamente e a segunda através de uma 
sequência de pontos intermédios. 

Tudo o que se refere aos atributos das imagens 
tais como luminosidade, movimento, tempo de expo- 
sição, etc. é dado através dos parâmetros das subro- 
tinas. 

A função de transformação da imagem para a 
sua colocação no ecran em coordenadas absolutas 
pode ser feita pelo «software» desde que o utilizador o 
especifique nos parâmetros correspondentes, Isso per- 
mitir-lhe-ã enviar os resultados nas dimensões pró- 
prias do problema sem se preocupar com as dimen- 
sões e coordenadas do ecran. 

Tem-se portanto (1): 


1.º tipo 
Organização da «display-file= 
OPENDI 


2.º tipo 
Construção de imagens 
PLOT 

CLOSEI 


As diferentes subrotinas permitem realizar as se- 
guintes funções: 


OPENDI — indicar o início de saídas gráficas e 
receber do utilizador indicações so- 
bre transformações das imagens a 
representar quanto a translações e 
homotetias das mesmas; 


OPENIM — indicar o início de uma imagem; 


CLOSEI — além da indicação de fim de imagem 
esta subrotina contém as indicações 
sobre o tempo de representação, so- 
bre a fotografia e sobre se é a úl- 
tima imagem a representar; 


PLOT — jndicar o novo ponto a representar 
e receber indicações sobre a inten- 
sidade luminosa da sua representa- 
ção; 

LINE — além da indicação das variáveis in- 


dexadas que contêm as diferentes 
coordenadas dos pontos da linha a 


() Para a descrição completa das subrotinas e seus parâmetros consultar o manual do utilizador, 


TÉCNICA 437 


at9 


representar, contém informação so- 
bre a intensidade luminosa da sua 
representação, 


O utilizador deve chamar a subrotina OPENDI 
para preparar a construção da «<Display-file» pelo 
computador e construir cada imagem por sucessivas 
chamadas das subrotinas PLOT e LINE. As diferen- 
tes chamadas de PLOT e LINE com vista a cons- 
truir uma imagem devem ser enquadradas entre a 
chamada da subrotina OPENIM que indica o início 
duma imagem e a chamada da subrotina CLOSEI 
que indica o seu fim, 


O utilizador pode construir tantas imagens quan- 
tas quiser, tendo apenas limitações referentes aos 
tempos de exposição que impõe e à dimensão de 
banda utilizada para registar a informação para o 
G. D. 


A sucessão típica de chamadas para o terminal 
gráfico é portanto dentro do programa normal do 
utilizador a seguinte: 


CALL OPENDI — abertura do G. D. 
CALL OPENIM — abertura da 1.º ima- 
gem 


Construção esmo 
da imagem — — 
com auxílio — — 
de PLOT e | CALL PLOT 


LINE — 
CALL LINE 
CALL CLOSEI — fecho da 1.º imagem. 
CALL OPENIM — abertura da enézi- 
ma imagem 
CALL PLOT 
Construção came 
da imagem = 
com auxílio 
de PLOT e | —— 
LINE CALL LINE 
CALL CLOSEI — fecho da enézima 
imagem e do G. D. 
280 


Entre a chamada de quaisquer 2 subrotinas do 
PGBD podem existir quaisquer outras instruções 
FORTRAN do programa de aplicação. 


4 — Conclusão 


A possibilidade de definir o tempo de exposição 
de uma imagem no ecran e de fazer variar o valor 
desse parâmetro permite a construção de desenho 
animado a partir de instantâneos no tempo. Embora 
ainda não tenha sido convenientemente explorada essa 
possibilidade tem interesse em problemas cuja solução 
não é calculada de forma contínua (ainda que dis- 
cretizada em pequenos intervalos de tempo) como é 
o caso das aplicações aqui apresentadas. 

Além disso para os utilizadores com vocação 
artística esta poderá ser uma via de realização de 
desenho animado. 

Esta facilidade estã pouco implementada noutros 
sistemas mas pode permitir interessantes explorações 
não só de tipo artístico como didático. 


IV — APLICAÇÃO DO PGBD A UM PROBLEMA DE 
MECÂNICA DE FLUIDOS 


1 — Primeira utilização 


A fim de se testar a qualidade dos filmes que 
poderiam ser obtidos através do sistema descrito an- 
teriormente, escolheu-se o programa SMAC [11], que 
resolve problemas de mecânica de fluidos em domi- 
nios bidimensionais, com superfície livre. 

Duas das características do método toraam-no 
particularmente adequado à feitura de filmes; a uti- 
lização de partículas sem massa, para definir a super- 
fície livre, possibilita a visualização do escoamento; 
e o facto de se tratar de um método explícito, gas- 
tando pouco tempo de cálculo entre a geração de 
duas imagens (em virtude de o número de imagens 
ter que ser elevado, tem de usar-se um passo no tem- 
po pequeno, não existindo portanto problemas de ins- 
tabilidade). 

Antes de se fazerem testes com escoamento de 
fluidos, fizeram-se alguns testes preliminares para a 
geração de imagens fixas, nomeadamente circunferên- 
cias. Foi necessário nessa altura fazer ligeiras alte- 
rações à subrotina DISP do PGED, usando-se buffers 
de tamanho variável e contendo um número inteiro 
de imagens, 

A decisão de se usar o programa SMAKC para 
estes testes preliminares em vez de usar outros pro- 
gramas do CTAMPUL resultou de se disporem já 
de resultados gráficos obtidos numa unidade de mi- 
crofilmes de Los Alamos e de ser assim mais fácil 
fazer comparações. 

O teste que se executou foi o «BROKEN DAM 
NO SLIP WALLS» (rebentamento de uma barragem, 
fluido viscoso) em que se observa o comportamento 
de uma coluna de liquido, contida inicialmente entre 
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duas paredes, a partir do instante em que uma delas 
«desaparece». Supõe-se que na parede que resta e 
no solo a velocidade do fluido é nula. 

Usaram-se 98 partículas para visualizar o escoa- 
mento (embora se tenham usado 1056 nos cálculos), 
com buffers de 588 bytes, contendo cada buffer duas 
imagens iguais. Contando com o gap de MB bytes, 
obtivemos portanto um ritmo de aproximadamente 


30 000 bytes/seg. ' 
26 imagens/seg,. E ) 
(5884-548) bytes/imagem | 


Foram realizados 2 filmes do mesmo problema, 
um com a duração de 5s que precisou de 20 min. de 
tempo de CPU no IBM 360/44 do CCUL e outro de 
dO s que precisou de 2 horas de CPU no mesmo com- 


Situação inicial do fluido 


ti 
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putador. O facto de o tempo de CPU não aumentar 
linearmente com o tempo de filme é devido ao passo 
no tempo ser mais pequeno no 2.º filme (sendo infe- 
rior em média o número de iterações necessário para 
o cálculo da pressão), 

Foram tiradas algumas fotografias ao osciloscô- 
pio, em fases sucessivas do colapso da massa de li- 
quido que dão uma ideia da evolução do processo, 
como a seguir se representa, 

Da observação dos filmes chegou-se à conclusão 
que este sistema de saida de resultados do computa- 
dor é uma excelente ferramenta de trabalho no do- 
mínio da Análise Numérica aplicada à Mecânica dos 
Fluidos, assumindo os resultados muito maior signi- 
ficação do que milhares de algarismos a representar 
campos de velocidade e posições de partículas, 
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C. D. U. 62.001.6(09) 
J. 3. DELGADO DOMINGOS 


CENTRO DE TERMODINAMICA APLICADA E MECÂNICA 
DE FLUIDOS DAS UNIVERSIDADES DE LISBOA 
(GC. T. A: M. F. UU. L) 
(Breve resenha sobre a sua criação e evolução) 


Técnica N.º 437 — LI — 12-1976, p. 133 a 142 


Apresenta-se uma breve resenha sobre o Centro de Ter- 
modinâmica Aplicada e Mecânica dos Fluidos das Universi- 
dades de Lisboa desde a sua criação até ao momento actual. 
Focam-se as vicissitudes mais salientes por que passou, a 
orientação adoptada pela sua Direcção e a inserção do Centro 
no panorama português da investigação. O trabalho realizado 
sintetiza-se na lista de publicações directamente originadas 
no Centro e seus colaboradores. 


C. D. U. 629 
J. J. DELGADO DOMINGOS 
RECURSOS NATURAIS, ECONOMIA E SOCIEDADE 
Técnica N.º 437 — LI — 12-1976, p. 143 a 157 


Na primeira parte do artigo examinam-se alguns aspectos 
da economia corrente e do planeamento e extraem-se algu- 
mas ilações para o caso nacional. Na segunda parte, apre- 
sentam-se alguns resultados da análise energética aplicados 
à Agricultura, à Indústria e aos Transportes e mostra-se como 
a crise do ambiente e a crise alimentar têm subjacente a 
crise da energia. Na terceira parte apontam-se algumas pre- 
visões acerca da «nova ordem econômica internacional». 


C. D. U. 620.9 
J. 3. DELGADO DOMINGOS 
O PROBLEMA ENERGÉTICO PORTUGUES 
Técnica Nº 437 — LI — 12-1976, p. 159 a 167 


Analisa-se sucintamente o problema energético português 
com dados existentes em 1975 e tendo em conta a situação 
geral do Pais quanto a capacidade tecnológica, financeira e 
de emprego. Os problemas energéticos são abordados tendo 
em conta a capacidade de refinação, os consumos nos Trans- 
portes, na Indústria e na Agricultura. Afloram-se possiveis 
alternativas energéticas, e alguns problemas subjacentes de 
natureza politica e social, As implicações da produção mono- 
polista de electricidade são referidas. A opção nuclear não 
se justifica no actual contexto. 
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C. D. U. 620.9:338.91 
J. J. DELGADO DOMINGOS 
ALTERNATIVAS ENERGÉTICAS 
Técnica N.º “37 — LI — 12.1976, p. 169 a 173 


Qualquer actividade humana tem um custo em energia. 
Analisar o dispêndio em energia de todas as actividades é 
objectivo da Análise Energética. Neste artigo consideram-se 
valores já conhecidos para as principais actividades e ex- 
traem-se as ilações para o caso de Portugal. Conclui-se com 
a enumeração de um plano possivel de diversificação energé- 
tica para o caso português, tendo em conta os problemas do 
emprego e do desenvolvimento económico. 


C. D. U. 621.47 
J, ], DELGADO DOMINGOS 
APROVEITAMENTO DA ENERGIA SOLAR EM PORTUGAL 
Técnica N.º 437 — LI — 12-1976, p. 175 a 181 


Faz-se uma sucinta descrição das aplicações da energia 
solar e analisam-se as possibilidades actuais do pais quanto 
ao lançamento em larga escala da utilização de energia solar. 
Apontam-se como primeira meta as utilizações da energia 
solar no aquecimento. Esboça-se um programa concreto de 
actuação. 


C. D. U. 62.001 
J. 3. DELGADO DOMINGOS 


PLANEAMENTO DA INVESTIGAÇÃO NO SECTOR 
DA ENGENHARIA MECÂNICA EM PORTUGAL 


Técnica N.º 437 — LI — 121976, p. 183 a 186 


Como os recursos para a investigação são sempre escassos 
analisa-se a estratégia da sua repartição pelos sectores fun- 
damentais a partir do seu efeito multiplicativo ou factor de 
impacto. Tal impacto tem em conta a estrutura actual da 
indústria mecânica do país e os factores essenciais que têm 
motivado o seu desenvolvimento. Finalmente mostra-se como 
a investigação em Termodinâmica Aplicada representa um 
sector chave para um desenvolvimento moderno a caminho 
da progressiva autonomização, 
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C. D. U. 620.9:338.81 
]. 1 DELGADO DOMINGOS 
ENERGY OPTIONS 


Tecnica Nº 437 — LI— 121976, po 169 a 173 

Every human activity has an energy cost. To establish 
those costs is the object of Energy Analisys. This paper 
considers values already known for energy costs in the main 
human activities and points the implications for Portugal 
Some points for a viable program of energy diversification 
in Portugal are refered, taking into account the problems of 
employment and economic development, 


C. D. U. 621,47 
J. J. DELGADO DOMINGOS 
SOLAR ENERGY UTILIZATION IN PORTUGAL 


Técnica N.º 437 — LI — 124976, p. 175 a 181 

à brief reviow of solar energy utilization is done, in 
connection with present capabilities of implementation on 
the portuguese market. It |s ahown that a first goal to 
achieve |s the use of solar energy in heating applications. 
A plan of action in Portugal is sketeched. 


CG. D. U. 62.00] 
J. 1. DELGADO DOMINGOS 


PLANNING RESEARCH ON MECHANICAL ENGINEERING 
IN PORTUGAL 
Técnica N.º 437 — LI — 121976. p 183 a 185 

Because research funds are always scarce, its allocation 
must be based on impact analysis. This analysis takes into 
account the present structure and development of the mecha: 
nical industries In Portugal and the main factor behind its 
development. It is shown how research in Applied Thermody- 
namics is a key sector for modern development in a way 
to progressiva autonomy, 
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C. D. U. 62.001.6[09) 


). 31. DELGADO DOMINGOS 


A BRIEF SUMMARY OF THE EVOLUTION OF C.T.AM.FUL. 
Técnica Nº 437 — L| — 121976, p. 133 a 142 


A brief summary of the evolution of C.T.AM.PUL (Cen- 
ter for Applied Thermodynamics and Fluid Mechanics of the 
Universities of Lisbon) is presented together with the main 
difficulties which had and has to face, the main options of 
research and the reasons behind are refered, A list of publi- 
cations gives a synthesis of the main work done up to now 


C. D. U. 62.9 
1. 1. DELGADO DOMINGOS 
MATURAL RESOURCES. ECONOMY AND SOCIETY 
Tecnica Nº 47 — LI — 12-1975, p. 143 a 157 


The paper is divided in three parts. The first deals with 
some aspects of usual economy and planning and some 
warning are drawn for the Portuguese situation. The second 
part considers energy analysis in food production, industry 
and transportation and is shown how the environmental and 
food crisis have the energy crisis as the driving force. Third 
part, based on the previous data, traces shortly, a prediction 
for the «new international order in economys. 


Cc. D. U. 620,9 
1. 1, DELGADO DOMINGOS 
THE PORTUGUESE ENERGY PROBLEM 
Tecnica Nº 437 — Li — 121976, p. 158 a 167 


The Portuguese energy problem is briefly reviewed with 
data available in 1975, taking into account the technological 
capacity, financial resources and employment. The analysis 
considers refining capacity, Transportion, Industry and Agri- 
culture. Energy alternatives and some underlying political and 
social problems are refered. The implications of monopolistic 
production of electricity are mentioned. Nuclear power stations 
are shown to have no justification on present situation, 
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C. D, U. 620,93:536.662 
MÁRIO NERY R. NINA 


FENÔMENOS DE COMBUSTÃO COM IMPORTÂNCIA 
EM ENGENHARIA 


Técnica N.º 437 — LI — 12-1976, p. 187 a 190 


Portugal depende dos combustiveis fósseis para satisfazer 
mais de B0%%, das suas necessidades energéticas presentes. 

Um conhecimento mais aprofundado da Mecânica dos 
Fluidos e da Quimica dos Sistemas com Combustão será 
indispensável se se pretender melhorar o aproveitamento 
dos combustiveis e minorar os efeitos prejudiciais das emis- 
sões que lhe estão associados. 

& linha de Combustão do C.T.A.M.F.U.L. propõe-se con- 
tribuir para este esforço, desenvolvendo métodos de cálculo 
de câmaras de combustão e validando-os através de expe- 
rimentação laboratorial e de medidas em instalações indus- 
triais, divulgando os resultados obtidos e realizando cursos 
intensivos na sequência do que efectuou em Dezembro de 1975. 


Cc. D. U. 621.018 


J. ]. DELGADO DOMINGOS 
E. MANSO 
MARIO NERY R. NINA 


FORNALHA EXPERIMENTAL CALORIMÉTRICA 
Técnica N.º 437 — LI— 12-1976, p. 191 a 194 


Projectou-se e construiu-se uma fornalha de secção qua- 
drada (500x500 mm) com o objectivo de permitir obter valores 
experimentais para validação de métodos de cálculo de cã- 
maras de combustão e fazer ensaios de queimadores e de 
corrosão de alta temperatura. 

A estrutura da fornalha permite variações na geometria 
dos circuitos dos gases e da água de arrefecimento. 

Apresentam-se resultados de ensaios de um queimador do 
tipo cone rotativo, onde se mostra a influência da carga e do 
excesso de ar na distribuição do fluxo de calor. 
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C. D. U. 536.24:620.197:697.8 
J. M. BRISSON LOPES 


TRANSFERENCIA DE CALOR EM REGIME ESTACIONARIO 
EM CHAMINES (ANALISE E RESOLUÇÃO NUMÉRICA) 


Técnica N. 437 — L| — 124576, p. 195 a 205 


O presente trabalho surgiu da necessidade do conheci- 
mento da temperatura e velocidade dos gases à saida duma 
chaminê conhecidos os seus valores na base. Tais valores 
são importantes para a previsão da dispersão de poluentes 
emitidos pelas chaminés para a atmosfera. Analisa-se a trans- 
ferência de calor entre os gases e o exterior das chaminés. 
Calculam-se os coeficientes de transmissão superficial de 
calor e as propriedades do ar a 760 mmHg. Realizou-se assim 
um programa para cálculo da temperatura dos gases ao longo 
da chaminé, da sua velocidade na base e no topo e os pertis 
de temperatura nas paredes interior e exterior. Analisam-se 
seguidamente alguns casos reais, sendo tratado o problema 
da corrosão de chaminés metálicas. 


C. D. U. 628.5 


M. TERESA PONTES 


FACTORES DO DIMENSIONAMENTO DUMA CHAMINE 
EMISSORA DE POLUENTES 


Técnica N.º da? — LI — 124576, p. 207 a 214 


Neste trabalho descreve-se sucintamente o cálculo da 
altura duma chaminé emissora de poluentes e acentuam-se 
os pontos criticos: concentrações máximas de poluentes a 
admitir no solo, condições meteorológicas e sobre-elevação 
da pluma, 

Como exemplo da sensibilidade do modelo às condições 
de cálculo apresenta-se o dimensionamento da chaminé prin- 
cipal da Refinaria da Petrosul em Sines. 


C. D. U. 536.2:65.046.3:664.8.037:621.58.336.24 


1. M. BRISSON LOPES 


SIMULAÇÃO NUMÉRICA DA CONGELAÇÃO DE PEIXE 
POR CIRCULAÇÃO FORÇADA DE AR 


Técnica N.º 437 — LI — 124976, p. 215 a 222 


Resolve-se um problema bidimensional de congelação de 
peixe por meio dum programa de cálculo numérico que 
permite obter a solução dum sistema de equações às deriva- 
das parciais do tipo: 


u, ô"u, Su, 
=Ta— Wig—— + E 
Bt gx2 By 
com condições fronteiras do tipo 
ou 
jB yu) Tê=-=o 


0 


pelo método ADI. 

Nesta resolução toma-se em conta que a congelação não 
se verifica a uma dada temperatura, mas sim numa determina- 
da gama de temperaturas e que é desigualmente distribuida, 
conduzindo a um problema com propriedades fisicas variando 
discretamente. Apresentam-se também gráficos com a dis- 
tribuição das isotérmicas que evidenciam o avanço da frente 
de congelação ao longo do tempo. 
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C. D, U. 596.24:620.197:607.8 
1. M. BRISSON LOPES 


STEADY STATE HEAT TRANSFER IN CHIMNEYS 
(ANALYSIS AND NUMERICAL TREATMENT) 


Técnica Nº d7— LI — 121976. p. 195 a 205 


The prevision of the dispersion of airborne chimney or 
stack effluents requires that the temperature and the speed 
of the gases at the top of the stack should be known. This 
paper presents a way of computing those values given the 
temperature at the bottom, the mass flow rate of the gases 
and the caracteristic stack dimensions. An analysis of the 
heat transfer between the flue gases and the ouiside of the 
stack is made. The heat transfer coeficients and the properties 
of eir at 760 mm of Hg are celculsted. A computer pro- 
gram was developped. This program computes the longitudinal 
temperature profiles of gases and of the inside and the 
outside stack surfaces, and the gas speed at the top and 
bottom of stacks. The last part of this paper analyses a few 
real cases and presents an example of prevision of minimum 
inside stack surface temperature, am important parameter for 
preventing acid corrosion in metalic stacks. 


C. D. U. 628.5 
M. TERESA PONTES 
DESIGN FACTORS OF A POLLUTANT EMITTING STACK 
Técnica N.º 497 — LI — 12-1976, p. 207 a 214 


We briefly describe the calculation of the height of a 
pollutant emission stack. The critical points are stressed: 
maximum ground pollutant concentration, meteorological con- 
ditions and plume rise. 

às an example of the model's sensivity to the Input 
variables, we present the calculation of the Petrosul Refinery's 
main stack, 


C. D,. U, 536.2:66.045.3:864.8.037:621.58:336.24 


J. M. BRISSON LOPES 
NUMERICAL SIMULATION OF FISH FREEZING BY CORCED 
AIR CONVECTION 
Técnica Nº 437 — LI — 12-1976, p, 215 a 222 
This paper describes the numerical solution of a fish 


freezing process by means of a computer program which 
solves sets of coupled equations of the kind: 


ou, Su, a=U 
— = Eq ==— ga Eq + Ê.. 
6t Ex2 sy: 
with bondary conditions obeying to: 
SL, | 
E— yu). TES 
an 


applying en ADI method. 

It was assumed that freezing occurs in an interval of 
temperature, the rate of freezing varying throughout this inter- 
val of temperatures, as it would happen when freezing sea 
water. The result is a problem with discrete variation of 
physical properties which is solved for an exemple of fish 
freezing by forced air convection. 

The results are also presented in graphics showing how 
the freezing front and the isothermal lines progress throughout 
the freezing process. 
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C. D. U. 620.93:536.662 
MARIO NERY BR. NINA 


COMBUSTION PHENOMENA AND THEIR APPLICATIONS 
IN ENGINEERING 


Técnica N.º 437 — LI — 12-1976, p. 187 a 190 


Over 80% of the energy used in Portugal today is derived 
from fossil fuels. 

A better understanding of the Fluid Mechanics and the 
Chemistry of Reacting Systems will be required in order to 
improve the efficiency of combustion processes and simulta- 
neously decrease the amount of harmful effluents. 

The Combustion Group of CT.AM.FU.L. is concentrating 
efforts in the development of better furnace prediction methods, 
validated through experimental data obtained both in the 
laboratory and im industrial instalations, and offering short 
specialised courses like the one that took place in Decem- 
ber 1975. 


C. D. U. 621.018 


!. 3). DELGADO DOMINGOS 
E. MANSO 
MARIO NERY R. NINA 


à CALORIMETRIC EXPERIMENTAL FURNACE 
Técnica N.º 43! — LI — 12-197%6, p. 19] a 194 


A 500x500 milimeter square section furnace was designed 
and built to provide data for the testing of prediction proce- 
dures applied to furnace design and to enable experiments 
to be carried on high temperature corrosion and burner design, 

The furnace structure allows the geometry of the gas 
passage and the cooling circuit to be assembled in several 
ways, 

Some preliminary test results are presented and show 
the effect of excess of air and load on local heat flux using 
a rotating cup type of burner. 
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C. D. U. 621.47:621.383 
LEOPOLDO J. M. GUIMARÃES 


AS CÉLULAS FOTOVOLTAICAS NO APROVEITAMENTO 
DA ENERGIA SOLAR 


Técnica N.º 437 — L| — 121976, p. 223 a 225 


A energia solar representa uma solução que se considera 
possivel para resolver em parte o grave problema de encon- 
trar fontes de energia que possam substituir em termos de 
futuro as disponibilidades energéticas clássicas que num 
futuro mais ou menos próximo se encontrarão esgotadas. 

Neste artigo tecem-se algumas considerações acerca do 
aproveitamento da energia solar, usando células solares as 
quais convertem a radiação solar directamente em energia 
eléctrica com uma razoável eficiência. 


C. D, U. 697.97:662.997 
LF. GC: RORIZ 


ELABORAÇÃO DE TABELAS DE TEMPERATURA SOL-AR 
PARA PAREDES, PARA A REGIÃO DE PORTUGAL 
CONTINENTAL 


Técnica Nº 437 — LI — 12-1976, p. 227 a 235 


De todos os métodos utilizados para o cálculo de cargas 
térmicas em Ar Condicionado, o da Temperatura Sol-Àr é o 
mais preciso e é largamente utilizado. Por estes motivos e 
dado que as tabelas existentes não se referem à situação 
climatológica de Portugal, elaborou o autor, tabelas de tem- 
peratura Sol-Ar próprias para es caracteristicas do nosso pais. 


OC D. U. 517.93:518.1 


J. 1], DELGADO DOMINGOS 
L. F. C. RORIZ 


APLICAÇÃO DO METODO EXPLICITO -IMPLICITO 

EM SISTEMAS DE EQUAÇÕES PARABÓLICAS, 

NÃO LINEARES, EM DOMINIOS DE GEOMETRIA 
ARBITRARIA 


Técnica N.º 437 — LI — 121976, p. 237 a 28 


Define-se o método explicito-implicito para a solução 
numérica de sistemas de equações às derivadas parciais para- 
bólicas não lineares em domínios de geometria arbitrária e 
apresentam-se alguns resultados de aplicação. As vantagens 
relativas ao método dos elementos finitos e aos de diferenças 
finitas implicitas são referidas quanto aos aspctos memória 
— tempo de cálculo. 
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C. D. U. 518.1 
HEITOR LOBATO GIRÃO PINA 
O METODO DOS ELEMENTOS FINITOS 
Técnica N.º 437 — LI — 121976, p. 241 a 248 


Neste artigo faz-se uma introdução ao método dos ele- 
mentos finitos, descrevendo as suas principais caracteristicas 
que o tornam distinto dos outros métodos de solução de 
equações diferenciais. 


C. D. U. 621.154 
A. F. DE O. FALCÃO 


O METODO DAS SINGULARIDADES NO ESTUDO 
TRI-DIMENSIONAL DO ESCOAMENTO 
EM TURBOMAQUINAS AXIAIS 


Tésnica Nº 437 — L| — 12-1976, p. 249 a 256 


Apresentam-se os fundamentos do método das singulari- 
dades para representação do escoamento em turbomáquinas 
axiais. Em particular, descreve-se a aplicação dos métodos 
da linha e da superficie sustentadoras, para estudo do escoa- 
mento em cascatas rectilineas de pás com geometria tri-di- 
mensional. 


D. D. U. 621-154 


A. 3. DE ALMEIDA SARMENTO 
|. A. FALCÃO DE CAMPOS 


ESCOAMENTO POTENCIAL BIDIMENSIONAL 
NUMA CASCATA DE PERFIS RECTOS 
— MÉTODO DA TRANSFORMAÇÃO CONFORME 


Técnica N.º 437 — L| — 12-1976, p. 257 a 265 


Apresenta-se a análise do escoamento potencial bidimen- 
sional numa cascata inclinada de perfis rectos a partir da 
transformação conforme do escoamento em torno do circulo, 
Obtêm-se as expressões da velocidade, ângulo de ataque e 
forças que actuam sobre os perfis. Apresentam-se e discutem- 
-se os resultados numéricos obtidos para a distribuição de 
pressão na região da cascata, ângulo de deflexão, coeficiente 
de sustentação e coeficiente de Interferência. 
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C. D. U. 518.1 
HEITOR LOBATO GIRÃO PINA 
THE FINITE ELEMENT METHOD 
Técnica Nº 437 — LI— 121976, p. 241 a 248 


This paper presents the finite element method in an 
introductory way and describes its principal features which 
distinguishes it from other methods for solving differential 
equations. 


C. D. U. 621.154 
A. F. DE O. FALCÃO 


FLOW ANALYSIS IN AXIAL TURBOMACHINES 
By SINGULARITY METHODS 


Técnica N.º 437 — LI — 12-1976, p. 249 a 256 


The use of singularity methods for the analysis of the 
flow in axial turbomachines is outlined. In particular, lifting- 
“line and lifting-surface theories are employed to study the 
flow about a three-dimensional straight cascade of twisted 
blades. 


D. D. U. 621-154 


A. J. DE ALMEIDA SARMENTO 
J. A. FALCÃO DE CAMPOS 


TWO-DIMENSIONAL POTENTIAL FLOW THROUGH 
A CASCADE OF STRAIGHT BLADES 
— METHOD OF THE CONFORMAL TRANSFORMATION 


Técnica N.º 437 — L| — 12-1976, p. 257 a 265 


The theoretical analysis of the two-dimensional potential 
flow about a staggered cascade of straight blades is presented 
based on the conformal transformation of the flow about a 
circle. Expressions are derived for the velocity, angle of 
attack and forces acting on the aerofoils. Numerical results 
for the pressure distribution about the cascade, deflection 
angle, lift coefficient and interference coefficient are also 
presented and discussed. 


C. D. U. 621.47:621.383 
LEOPOLDO J. M. GUIMARÃES 


PHOTOVOLTAIC CELLS IN DIRECT SOLAR ENERGY 
CONVERSION 


Técnica N. 437 — LI — 12-1976, p. 223 a 225 


Solar energy represents one possible solution to the 
major problem we face today: how to find energy resources 
to replace the non-renewable energy reserves we have been 
using and wasting up to now. 

In this paper some considerations are made on solar 
cells which are photovoltaic cells that convert sunlight directly 
into electrical energy with a good efficiency. 


C. D. U. 697.97:662.997 
L. F. C. RORIZ 


SOL-AIR TEMPERATURE TABLES OF WALLS FOR PORTUGAL 
(CONTINENTAL) 


Técnica N.º 437 — LI — 12-1976, p. 227 a 235 


The Sol-Air temperature method is an accurate and widely 
used method in àAir-Conditoning engineering, for the cal- 
culation of heat gains in a room, This work presents tabulated 
values that are to be used in Portugal, as evaluated tables of 
Sol-Air temperature for this weather conditions are not 
available. 


C. D. U. 517.93:518.1 


3. 31. DELGADO DOMINGOS 
L. F. C. RORIZ 


IMPLICIT-EXPLICIT METHOD APPLIED TO NON LINEAR 
PARABOLIC DIFFERENCE EQUATIONS IN A ARBITRARY 
GEOMETRY 


Técnica N.º 437 — L| — 12-1976, p. 237 a 239 


The explicit-implicit method for the numerical solution of 
systems of non linear parabolic partial differential equations 
in arbitrary geometries is defined and some results of applica- 
tion are shown. The relative advantages in relation to the 
finite element methods and implicit finite difference methods 
are refered in what concerns memory-computing time. 
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C. D. U. 62.002.1 
L. F. OC. RORIZ 


PROGRAMA PARA RESOLUÇÃO DUMA REDE 
DE ACTIVIDADES PELO METODO DO CAMINHO CRÍTICO 


Técnica N.º 47 — L| — 121975, p. 267 a 20 


Embora existam já programas para resolução dum diagra- 
ma CPM-tempo por computador, a sua utilização além de 
onerosa não era didáctica. Por estas razões foi desenvolvido 
no CTAMFUL, para apoio da Divisão de Termodinâmica Apli- 
cada, um programa para a resolução deste tipo de problemas. 
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LI — Dezembro 1976 


J, 1 DELGADO DOMINGOS 
J. ). AMARANTE DOS SANTOS 
M. QUIRINO 
A. NEVES DE CARVALHO 


PROJECTO, CONSTRUÇÃO E UTILIZAÇÃO 
DE UM «DISPLAY. ÓPTICO «OFF-LINE» 


Técnica Nº 47 — L| — 12-1976, p. 271 a 282 


A necessidade de representar de modo cómodo o enorme 
volume de dados gerados na solução numérica das equações 
às derivadas parciais que regem o escoamento de fluidos e o 
elevado custo do equipamento existente no mercado levou 
ao projecto e construção de um «display» óptico e ao desen- 
volvimento do correspondente «software» de operação. É este 
o trabalho que se descreve no artigo. O equipamento reali- 
zado, centrado sobre problemas de aplicação, além do seu 
baixo custo, permitiu evidenciar as amplas possibilidades que 
tal «display» oferece em inúmeros outros dominios. 
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DESIGN, MANUFACTURE AND USE OF AN OFF-LINE 
OPTICAL DISPLAY 


Técnica N.º 437 — LI — 12-1976, p. 271 a 282 


The need for display representation processing numerical 
data generated on the numerical solution of partial differential 
equations which govern the flow of viscous fluids, and the 
very expensive equipment available on the market for visual 
display motivated the development of the prototype whose 
hardware and software are described. 

The work done stressed the importance of project deve- 
lopment centered on applications and showed the cheapness 
and wide range of applicability of the equipment built. 
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C. D. U. 62.002.1 
L. F. C. RORIZ 
PROGRAM FOR SOLVING CPM-TIME 
Técnica N.º 437 — LI — 121976, p. 267 a 270 


None of the existing programs for solving CPM-time 
problems is suitable for teaching purposes, to the author's 
best knowledge as they are all available in machine language 
only. Therefore, a new CPM-time solving program was deve- 
loped at CTAMFUL, with teaching purposes in mind, and has 
already been used by the Applied Thermodynamics Department 
students. 
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ERRATA 


SIMULAÇÃO NUMÉRICA DA CONGELAÇÃO DE PEIXE 
POR CIRCULAÇÃO FORÇADA DE AR 


ONDE SE LÊ 


pág. 217 — 1.º coluna 


10º linha | 


| pãg. 219 — 1.º coluna 
7." linha 


10." linha 


19.º linha 


DEVE LER-SE 
ôu, 
Pe 18) M$ = 


= 42244812 erg gi K-1 
= (0,487 kcal h2 ma K-1 = 56 597,2 erg s1 cmi Ka 


480 X 252 X 


+ ii 


DR 


CÁLCULO 
DE 


POÓRTICOS 


Método de cálculo simples e rápido, entrando 
em linha de conta com o deslocamento dos nós 
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